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FOREWORD 

The  studies  reported  here  were  conducted  largely 
under  U.  S.  Air  Force  Contract  AF49(638)-817 , 
monitored  by  the  Air  Force  Office  of  Scientific 
Research,  of  the  Office  of  Aerospace  Sciences. 

The  analytical  work  was  incomplete  at  the  end  of 
that  contract.  The  final  formulation  of  the  com¬ 
bustion  model  was  completed  with  funds  from  a 
Rocketdyne  Liquid  Propellant  IR&D  program  and  from 
Contract  NAS  8-19  (J-2  program).  The  analytical 
formulation  was  put  into  a  form  useful  for  pre¬ 
dicting  combustion  characteristics  of  rocket 
engines  of  interest  in  those  programs. 
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ABSTRACT 


Simultaneous  equations  describing  rocket  propellant 
injection,  atomization,  mixing,  vaporization,  and 
combustion  are  formulated  for  a  cylindrical  liquid 
oxygen  jet  surrounded  by  an  annular  gaseous  hydrogen 
stream.  The  formulation  is  based  on  division  of  the 
combustion  chamber  into  a  nonburning  region  near  the 
injector  and  a  combustion  region  further  downstream. 
Processes  in  the  nonburning  region  are  required  to 
satisfy  certain  combustibility  criteria  before  com¬ 
bustion  is  possible.  The  propellents  then  pass 
through  a  plane  flame  front  into  the  combustion 
region.  A  simplified  treatment  of  turbulent  mixing 
betweenunlike  reacting  gas  streams,  one  of  which 
contains  a  liquid  oxygen  spray,  characterizes  the 
combustion  region  model.  Solutions  of  the  system 
of  equations,  obtained  with  a  digital  computer  pro¬ 
gram,  are  detailed,  discussed,  and  compared  favor¬ 
ably  with  experimental  information  from  Part  I. 
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SUMMARY 


A  complete  analysis  of  the  steady-state  rocket  combustion  processes 
attending  coaxial  jet  injection  of  liquid  oxygen  and  gaseous  hydrogen 
is  presented.  Programmed  for  IBM  709^  ID  Ml  solution,  a  model  is  formu¬ 
lated  that  consists  of  simultaneous  equations  describing  LOX-jet  atomiza¬ 
tion,  gas-spray  mixing,  spray  vaporization,  spray  ballistic  behavior, 
gas  stream  expansion  and  deceleration,  gas  stream  combustibility,  pro¬ 
pellant  combustion,  turbulent  mixing  between  parallel  gas  streams,  and 
gas  dynamics  of  the  subsonic  portions  of  the  combustor.  Required  model 
input  parameters  are  those  normally  known:  injector  face  area,  oxidizer 
and  fuel  injection  areas,  flowrates  and  temperatures,  chamber  pressure, 
and  combustion  chamber  geometry. 

There  are  two  discrete  parts  to  the  model.  The  first  is  adjacent  to 
the  injector  and  is  characterized  by  the  absence  of  propellant  combustion. 
A  portion  of  the  liquid  oxygen  jet  is  atomized  and  mixed  in  a  controlled 
manner  with  the  surrounding  hydrogen  stream;  a  cyl indrically  annular, 
layered  spray  structure  is  hypothesized  with  the  earliest  formed  spray 
occupying  the  outer  layers.  The  hydrogen  stream  is  decelerated  and  ex¬ 
panded  radially  by  diffusion  and  mixing  with  the  surrounding  (stagnant) 
gases  and  by  momentum  exchange  with  the  contained  liquid  oxygen  spray. 
Spray  vaporization  effects  a  gradually  increasing  oxygen  concentration 
in  the  hvdrogen  stream.  Recirculation  of  combustion  gases  from  the 
combustion  region  that  lies  downstream  is  presumed  not  to  occur;  support¬ 
ing  evidence  is  cited. 

The  gases  in  the  nonburning  region  are  considered  to  be  combustible  if 
their  hydrogen  concentration  is  lower  than  the  upper  flammability  limit 
and  if  their  velocity  is  lower  than  the  local  turbulent  flame  speed. 
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When  these  conditions  are  satisfied,  flame  spreading  is  rapid  enough 
that  the  nonburning  region  is  terminated  by  a  plane  flame  front  standing 
in  the  spray- laden  gases. 

The  second  part  of  the  model  is  a  combustion  region.  A  portion  of  the 
gases  in  the  nonburning  spray  zone  is  assumed  to  be  diverted  around  the 
flame  front  and  form  a  fuel-rich  gaseous  mantle  around  a  central,  well- 
mixed,  one-dimensional  spray  combustion  zone.  The  progress  of  propellant 
combustion  can  be  limited  by  these  processes:  (l)  atomization  of  the 

residual  liquid  oxygen  jet  that  penetrates  into  the  combustion  zone,  (2) 
vaporization  of  liquid  oxygen  droplets,  and  (3)  turbulent  mixing  of  the 
gaseous  mantle  into  the  spray  combustion  zone. 

Few  data  exist  concerning  some  of  the  processes,  so  the  model  formulation 
includes  a  number  of  arbitrary  constants  and  assumptions.  Selected  values 
for  these  were  modified  until  a  combination  was  found  that  would  permit 
the  model  predictions  to  reproduce  experimental  information  from  two 
particular  transparent  model  rocket  motor  tests.  The  results  from  those 
best-fit  model  calculations  are  presented  in  enough  detail  to  characterize 
prominant  features  of  the  model  predictions. 

Although  the  best-fit  cases  agree  well  with  those  two  tests  which  had 
given  the  most  complete  experimental  data,  parametric  variations  of  model 
chamber  pressure,  mixture  ratio,  and  propellant  injection  temperatures 
gave  some  predictions  that  did  not  agree  particularly  well  with  observed 
trends  in  other  tests  which  had  given  incomplete  experimental  data. 

These  comf prisons  are  discussed  in  terms  of  the  model's  validity  and 
additional  definitive  measurements  that  might  have  been  made  in  the  tests. 
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INTRODUCTION 


This  report  is  the  second  in  a  series  presenting  the  results  of  a  combined 
experimental  and  analytical  study  of  the  processes  that  transpire  during 
steady-state  rocket  combustion  with  liquid  oxygen  and  gaseous  hydrogen 
propellants.  The  first  report  in  the  series  detailed  the  experimental 
investigation  with  three  injector  types  and  included  a  qualitative  dis¬ 
cussion  of  the  combustion  processes  found  to  be  important  with  each  type 

(Hrf.  l).  As  noted  there,  specific  experimental  information  was  sought 
to  aid  in  formulating  physically  realistic  analytical  models  of  the  bi- 
propellant  combustion  processes.  The  analysis  that  has  been  based  on  the 
experimental  information  concerning  coaxial  jet  injection,  in  which  each 
injection  element  consists  of  u  cylindrical  liquid  oxygen  jet  surrounded 
by  an  annular  gaseous  hydrogen  stream,  is  the  subject  of  this  report. 

The  experiments,  conducted  in  a  transparent-walled  model  rocket  motor, 
revealed  that  the  liquid  oxygen  jets  penetrated  about  3  to  4-1/2  inches 
downstream  of  the  injector  where  they  rather  abruptly  ceased  to  be  visible 
in  high-speed  motion  pictures.  At  about  the  same  position,  streak  photo¬ 
graphy  indicated  an  abrupt  increase  in  the  velocity  of  luminous  combustion 
elements  presumed  to  be  representative  of  combustion  product  gases.  A 
short  distance  further  downstream,  melting  of  the  pyrex  combustor  walls 
began.  These  phenomena  were  interpreted  to  result  from  a  flame-front 
standing  at  a  particular  position  in  the  oxygen-spray-laden  gas  flow 
surrounding  each  liquid  oxygen  jet.  A  derived  schematic  representation 
of  the  atomization,  jet-spreading,  and  combustion  processes  for  a  single 
coaxial  jet  injection  element  is  shown  in  Fig.  1  . 

The  analyses  presented  and  discussed  are  almost  exclusively  based  on  the 
concepts  shown  in  Fig.  1  .  Immediately  downstream  of  the  injector,  the 
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initial  liquid  oxygen  jet  atomization,  spray  acceleration  and  evaporation, 
and  hydrogen  stream  spreading  and  deceleration  are  presumed  to  take  place 
in  the  absence  of  direct  combustion.  Upstream  recirculation  of  combustion 
gases  is  assumed  either  not  to  occur  or  to  be  so  weak  as  to  have  no  appre¬ 
ciable  affect  on  the  important  processes  in  this  region,  A  resultant 
natural  division  of  the  model  formulation  for  this  injection  type  into  a 
distinct  nonburning  injection  region  and  a  subsequent,  downstream  com¬ 
bustion  region  is  used. 

A  complete,  integrated  combustion  model  for  the  entire  combustion  chamber, 
from  the  injector  face  to  the  nozzle  throat,  was  the  goal  pursued  in  this 
investigation.  At  many  points  in  the  model  formulation  it  became  neces¬ 
sary  to  select  one  approach  or  assumption  from  among  alternatives.  In 
some  instances,  review  of  the  experimental  data  suggested  the  most  real¬ 
istic  choice;  in  others,  peripheral  analyses  were  performed  to  guide  the 
development.  Occasionally,  a  choice  was  made  that  was  later  abandoned 
because  results  were  computed  which  could  not  be  supported  by  the  ex¬ 
perimental  evicence.  The  mo  -e  important  of  these  alternatives  are 
indicated . 
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MODEL  FORMULATION 


Injectors  for  liquid  propellant  rocket  combustors  are  usually  based  on  a 
single  fundamental  injection  concept  with  only  minor  modifications  or 
variations  to  provide  for  special  requirements,  such  as  film-coolant  *or 
chamber  walls  and  adequate  ignition  sources.  Thus,  a  high-thrust  rocket 
engine  injector  might  have  several , hundred  coaxial  jet  injection  elements, 
all  having  identical  geometric  characteristics  and  ostensibly  having 
identical  propellant  supply  and  flow  characteristics,  uniformly  spaced 
in  the  available  injector  face  area.  In  the  ensuing  formulation,  it  is 
assumed  that  the  processes  downstream  of  any  one  of  those  elements  may 
be  represented  by  a  stream  tube  which  is  characteristic  of  the  entire 
combustion  field.  The  model,  therefore,  is  based  on  division  of  the 
total  injection  area,  chamber  cross-sectional  area  along  the  chamber 
length  and  propellant  flowrates  by  the  number  of  elements  to  get  equiv¬ 
alent  values  per  element.  An  additional  simplification  in  this  vein  is 
the  assumption  of  cylindrical  symmetry  of  the  stream  tube  about  the 
element's  axial  centerline. 

The  model  formulation  consists  of  determining  the  important  processes  and 
developing  analytical  or  empirical  expressions  to  describe  them.  Since 
a  large  number  of  interrelated  phenomena  are  accounted  for,  the  system  of 
equations  and  expressions  must  be  solved  numerically.  An  IBM  7094  EDfM 
program,  written  in  FORTRAN  II  and  based  on  simple  stepwise  progression 
down  the  length  of  the  chamber,  has  been  employed;  conditions  at  a  plane 
x  distant  from  the  injector  were  used  to  calculate  conditions  at  a  plane 
(x  +  A  x)  distant  from  the  injector.  It  was  assumed  that  small  enough 
Ax  increments  could  be  used  to  avoid  mathematical  instability  and  accumu¬ 
lation  of  significant  errors  along  the  chamber  length. 
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NONBURNING  INJECTION  REGION 


Under  the  experimental  conditions  detailed  in  Ref.  1  ,  coaxial  jet  in¬ 
jection  of  gaseous  hydrogen  and  liquid  oxygen  was  characterized  by  a 
relatively  low-velocity  liquid  oxygen  stream  being  completely  surrounded 
by  a  much  higher  velocity  gaseous  hydrogen  flow.  At  injection-end 
chamber  pressures  of  **50  to  525  psia,  liquid  oxygen  injection  velocities 
ranged  from  about  50  to  100  ft/sec  (depending  upon  chamber  pressure  and 
mixture  ratio),  while  hydrogen  injection  velocities  ranged  from  about 
1000  to  3000  ft/sec  (depending  upon  chamber  pressure,  mixture  ratio,  and 
injection  temperature).  The  hydrogen  stream  formed  a  cylindrical  barrier, 
for  at  least  a  short  distance  from  the  injector,  which  insulated  the  con¬ 
tained  liquid  oxygen  stream  from  the  surroundings.  The  initial  liquid 
jet  atomization,  spray  acceleration,  and  spray  evaporation  unquestion¬ 
ably  took  place  in  the  absence  of  combustion.  This  section  is  concerned 
v  i  tli  the  model  formulation  for  that  nonburning  region. 

Liquid  Oxygen  Jet  Atomization 

Immediately  upon  leaving  the  end  of  its  confining  tube,  an  elemental 
cylindrical  plug  of  liquid  oxygen  finds  its  surface  exposed  to  tremendous 
shear  forces  from  the  surrounding  high-velocity  hydrogen  stream.  These 
forces  rapidly  induce  capillary  wave  motion  on  the  liquid  jet  surface. 
Mayer  (kef.  2  )  has  extended  an  analysis  of  capillary  waves  on  plane, 
deep  liquids  to  a  formulation  for  liquid  jet  atomization  processes.  He 
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derived  an  equation  relating  a  mean  droplet  diameter  produced  by  wave 
amplification  and  crest  shedding  to  the  gas  and  liquijd  stream  properties:* 


D 


mean 


=  (9)(l2)(l6)->  n  B 


TT  2 

D  U 
-  Hg  r 


/*>1/8 


2/3 


(1) 


where  B^  is  a  composite  numerical  factor  including  a  sheltering  parameter 
and  a  parameter  dependent  upon  the  crest  configuration  at  the  instant  of 
erosion  or  shedding.  On  conceptual  grounds,  Mayer  reasoned  that  should 
be  of  the  order  of  unity.  He  made  a  favorable  parametric  comparison  of 
his  theoretical  result  with  experimental  data  (Hef.  J  )  on  costream  and 
contrastream  atomization  of  cylindrical  molten  wax  jets  and  concluded 
with  an  empirical  evaluation  from  that  data  of  =  0,3.  This  was  further 
resolved  to  indicate  that  about  l/7  of  the  wave  amplitude  is  shed  as  a 
ligament  at  the  instant  of  erosion  and  that  about  70  percent  of  the  wave 
amplitude  is  sheltered  from  the  gas  stream. 


Mayer's  droplet  size  expression  has  been  used  in  the  model  formulation  to 
calculate  the  mean  LOX  droplet  diameter  being  produced  in  the  nth  Ax 
increment  from  the  injector:** 


r 


«  2719.2  Ba 


<V>1 


<v> 


(la) 


*A11  equations  have  been  modified  to  be  consistent  with  the  units  given 
in  the  nomenclature.  Values  of  the  required  properties  appear  in 
Appendix  A. 

**0r  of  the  nth  Bpray  element  for  instances  where  the  atomization 
process  ceases  for  a  time  and  is  then  resumed. 
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In  contrast  with  Mayer's  exposition  in  Ref.  2,  where  B^  was  evaluated 
from  the  entire  spray  resulting  from  complete  jet  atomization,  Eq.  la 
clearly  denotes  a  mean  drop  size  for  the  spray  produced  in  a  small  local 
region.  For  that  reason,  th  value  of  was  left  as  an  arbitrary  con¬ 
stant  in  the  model  formulation. 


Dickerson  and  Schuman  (Ref.  4)  have  shown  that  Mayer's  theory  can  be 
extended  to  provide  prediction  of  a  spray  production  rate.  Their  deriva¬ 
tion  was  in  terms  of  mass  atomization  rate  per  unit  liquid  surface  area: 


where  is  a  numerical  factor  related  to  in  Eq.  1. 

For  a  cylindrical  jet  segment  Ax  long,  the  local  spray  production  rate 
of  the  nth  spray  group  is  here  taken  as: 

(WAn)l_2  *  (6-859 7  I  I0*3)  CA  D{|4i  F(n) 

(2q) 

in  which  a  suppression  factor  F(n)  has  been  introduced  to  account  for  the 
combined  effects  of  finite  wave  buildup  time  and  finite  tube  wall  thickness 
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in  delaying  the  atomizalion  immediately  downstream  of  the  injector.  A 
parabolic  form  was  assumed: 


F(..) 


;  1  £  n  <  c2 


(3) 


Values  of  c^  and  were  usually  selected  to  suppress  the  atomization 
by  80  to  90  percent  at  the  injector  face  and  continue  suppression  for 
about  l/2  inch  from  the  injector. 

The  mass  of  an  individual  droplet  was  given  by 


Wdn  7)  ^dn  ^dn 


so  the  number  flowrate  of  droplets  in  the  nth  spray  group  was: 


Nd„  -  <V,  /wd„i 

1  "fc 


(*) 


(5) 


According  to  Kq.  1  and  2,  local  spray  droplet  size  increases  and  local 
spray  production  rate  decreases  with  decreasing  relative  velocity  between 
the  liquid  jet  and  the  atomizing  gas  stream.  As  will  become  aj^arent 
later,  momentum  exchange  from  the  high-velocity  gas  stream  to  the  acceler¬ 
ating  spray  can  cause  fairly  rapid  deceleration  of  the  gases  which,  in 
turn,  is  followed  by  production  of  fewer,  but  larger,  droplets.  Very 
large  dro  plets  would  later  be  accelerated  so  slowly,  however,  that  their 
velocities  over  much  of  the  chamber  length  might  be  lower  than  the  mini¬ 
mum  velocities  reduced  from  the  experimental  streak  photographs.  There¬ 
fore,  in  the  model  formulation,  a  maximum  droplet  diameter  was  imposed 
on  the  nonburning  atomization  process;  the  maximum  size  was  selected  by 
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n  peripheral  analysis  described  in  Appendix  B.  Calculation  of  droplet 
sizes  larger  than  the  assumed  maximum  was  equivalent  to  cessation  of  the 
atomization  process;  generally,  atomization  would  not  be  resumed  until  a 
position  in  the  chamber  was  reached  where  propellant  combustion  resulted 
in  a  relatively  high  gas  velocity. 


Individual  Droplet  Behavior 

The  behavior  of  an  individual  droplet  of  liquid  oxygen  in  a  convective 
gas  stream  is  strongly  dependent  upon  the  gas  stream  properties  and  con¬ 
ditions.  Three  kinds  of  behavior  were  accounted  for  in  the  model  formu¬ 
lation,  viz.,  ballistics  (acceleration  or  deceleration),  heating  (droplet 
temperature  clianges),  and  vaporization.  Deviations  from  spherical  shape, 
ii  terference  or  reagglomerat ion  with  neighboring  droplets,  and  aerodynamic 
breakup  of  a  droplet  following  its  initial  formation  were  not  considered 
to  be  of  significance  to  the  nonburning  analysis. 

Droplet  ballistics  were  described  by  the  often-used  drag  equation  for  a 
sphere  with  a  low  mass  evaporation  rate: 


(^12>CDPgl)rlUr 
"dt  '*  0d  Dd 


(6) 


Defining  a  vaiiable  K 
*  n 


for  the  nth  droplet  group,* 


3  CDnl  *  1 

**  ^dnl  ^dnl 


(7) 


♦Where  gas  properties  are  n-subscripted,  those  gases  locally  associatjd 
with  the  nth  spray  group  are  meant.  Defer  to  the  formulation  beginning 
on  page  21  for  a  description  of  this  a-  ociation. 
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Eq.  6  was  integrated  to  obtain: 


Ud„2  -  °-5  <<«„♦  »Ud„l  -KnUg„l  + 


K 

n 


U«rnl  '  (Kn  +  1)  UHnl  +  4  K  U  T  (U 

gnl  n  dnl  n  irnl  '  i 


n  gnl  gnl 


(8») 


For  rapidly  decelerating  droplets,  Eq.  8a  may  have  a  negative  radicand, 

so  the  less  accurate  form,  previously  used  in  lief.  5»  was  employed  when 

U  .  <  U  .  . : 
gnl  dnl 


U .  0  =  U .  .  + 

dn2  dnl 


K  (U  _  -  U.  .)  (U  .  -  U.  J 
n  x  gnl  dnl7  N  gnl  dnl 


U 


(8b) 


dnl 


The  composite  form  of  the  droplet  drag  coefficient  that  accounts  for  drop¬ 
let  flattening  at  high  droplet  Reynolds  numbers  (Ref.  5)  was  used: 


CDnl 

-  2,1  (“'dnl)'0'8' 

, Re ,  .  <  80 
dnl 

CDnl 

'  °’271  ('"dnl)0'"17 

,80  £  Re  .  ,  <  104 
dnl 

CDnl 

_  O 

'H'dnl  2  104 

(9) 

Free-stream  gas  properties  were  used  in  evaluating  the  Reynolds  numbers 
for  drag  coefficient  calculations. 
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Changes  in  droplet  bulk  temperature  were  calculated  by  means  of  equations 
resulting  from  analysis  of  simultaneous  heat  transfer  to  and  mass  transfer 
from  an  evaporating  droplet  (Ref.  6): 


T  =  T  + 
dn2  *  ’ 


144  A  x 


dnl  ,  Pdm  D-dnl  "am  + 


ilnl 


kfnl  Znl  NuIInl  ^Tfnl  '  Tdnl^ 


3600  (e  nl  -  1) 


<1M>  MT  "vnl  P1  NuXnl  .  f  ■'! 

- * — rrs -  1oe_  — 


RT 


f  nl 


p  _  p 
1  vnl. 


(10) 


where 


(1W)(3600)(C  )_1  M,  Fj  N'u^,  x  p 


nl 


_ vnl _ 

11  Tfnl  kf.,l  Nullnl 


log 


^  *1  '  Fvnl. 


(U) 


Here  the  Reynolds,  Prandtl ,  Schmidt,  and  Nusselt  numbers  were  evaluated 
using  approximations  of  the  droplet  film  properties  based  on  an  arithmetic 
mean  temperature  between  the  droplet  and  gas  temperatures. 


For  low  evaporation  rates,  the  change  in  liquid  droplet  mass  across  an 
increment  is  given  (Ref.  7)  by: 


2tt  M 


0f)  ^gnl  *vnl  *dnl  ^UMnl 


A  x 


U  =  W 

dr>2  dnl 


Vtk  Mgnl  ^  dnl  +  1  dn2) 


P  -  X  P 

vnl  0Qgnl  1 

P  -  P 
1  vnl 


(12) 
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The  term  in  brackets  could  become  negative  if:  (l)  the  droplet  tempera¬ 
ture  (and  thus  vapor  pressure)  were  very  low,  (2)  the  chamber  pressure 
were  ver>  high,  or  (l)  the  oxygen  concentration  in  the  gas  stream  were 
high.  The  digital  computer  program  formulation  avoided  the  implied  oxygen 
condensation  by  not  permitting  w  >  v(jn|*  Modifications  of  Fq .  12  to 
account  for  high  mass  transfer  rates  (Ref.  7)  were  not  considered  to  be 
necessary  in  the  nonburning  region. 

Evaporation  from  the  liquid  oxygen  ,jet  was  approximated  bv  assuming  that 
its  vaporization  rate  per  unit  surface  area  was  equal  to  that  of  the 
last-formed  spray  group. 


Hydrogen  Stream  Expansion  and  Deceleration 

At  a  position  just  downstream  of  the  injector,  the  propellunt  streams  may 
be  idealized  as  two  concentric  cne-d imens ional  flows.  They  immediately 
interact  with  each  other  and  with  the  surrounding  chamber  gases,  however, 
to  form  an  exceedingly  complex  bipropellant  flow  system.  The  atomization 
process  described  earlier  is  a  portion  of  the  interaction.  Now,  the 
interactions  of  the  hydrogen  stream  with  the  exterior  surroundings  and 
with  the  interior  liquid  oxygen  spray  will  be  considered.  Occasional 
reference  to  Fig.  2  may  be  helpful  in  understanding  the  hypotheses 
employed. 

The  outer  periphery  of  the  annular  gaseous  hydrogen  stream  exerts  viscous 
shear  forces  on  the  surrounding  chamber  gases  setting  them  in  motion 
alongside  the  high-velocity  hydrogen  stream.  Injection  stream  turbulence 
and  shear-induced  turbulence  lead  to  fairly  rapid  gas  mixing  around  the 
hydrogen  stream.  In  that  mixing  zone,  momentum  exchanger  decelerates  the 
hydrogen  and  accelerates  the  chamber  gases.  The  mixing  zone  becomes 
thicker,  and  the  hydrogen  velocity  defect  associated  with  it  penetrates 
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further  into  the  hydrogen  stream  with  increasing  distance  from  the  injec¬ 
tor.  Except  for  the  liquid  jet  spray  core  (which,  as  discussed  later, 
makes  the  problem  considerably  more  complex  by  transferring  most  of  the 
deplet  ion  in  hydrogen  momentum  to  the  interior),  the  hydrogen  stream 
behavior  may  be  related  to  the  classical  gas  dynamic  problems  of  the 
free  (unconfined)  turbulent  circular  jet  and  the  bounded  (confined) 
turbulent  circular  jet. 

With  botli  free  and  bounded  gas  jets  having  essentially  rectangular  in¬ 
jection  velocity  profiles  (fully  developed  turbulent  flow),  the  thicken¬ 
ing  mixing  zone  involves  more  and  more  of  the  jet  stream  until,  at  some 
finite  distance  from  the  injector,  it  has  grjwn  inward  to  the  jet  center¬ 
line.  Denoting  that  downstream  position  as  x  ,  it  is  noted  that  a  dimin¬ 
ishing  core  of  unmixed  jet  stream  penetrates  at  injection  velocity  to  x^ 
and  there  vanishes.  The  gus  dynamics  downstream  of  x  huve  been  success¬ 
fully  described  theoretically  for  the  free  jet  (lief.  8)  und  studied 
extensively  for  bounded  jets  (e.g.,  Ref.  9).  The  bounded  jet  problem 
lias  not  yet  yielded  to  complete  analytical  description.  Some  results 
for  the  free  jet  are  (lief.  8):  (l)  the  maximum  jet  velocity  occurs  on 

the  jet  centerline  and  decuys  inversely  with  axial  distance  from  the  jet 
source,  (2)  the  jet  width  increases  linearly  with  that  distance,  (3)  the 
radial  profile  of  axial  velocity  has  a  characteristic  shape  similar  to 
a  Gaussian  distribution,  and  (4)  the  total  volume  of  fluid  in  motion 

increases  linearly  with  distance  from  x  .  The  bounded  jet  results 

P 

differ  from  the  free  jet  in  that:  (l)  the  velocity  at  large  distance 
downstream  from  x^  approaches  a  constant  one-dimensional  velocity  (pre¬ 
scribed  by  mass  continuity)  rather  than  approaching  zero,  (2)  momentum 
exchange  ut  the  jet  boundary  ingests  gases  which  are  not  freely  avail¬ 
able  from  the  surroundings  but  are  circulated  upstream  from  the  neighbor¬ 
hood  where  the  jet  fills  the  boundary  channel  (thus  forming  an  annular 
eddy),  and  (3)  the  total  mass  flowrate  is  limited  to  the  jet  injection 
f 1 owrate. 
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The  phenomena  upstream  of  x^  are  generally  well  known  from  experimental 
studies  but  have  not  confirmed  a  general  theoretical  treatment.  Even 
without  the  additional  complication  of  having  a  liquid  core,  it  is 
apparent  that  empirical  approximations  are  required  for  a  flow  model 
formulation  beginning  at  the  injector  face. 

Before  stating  the  approximations  used  for  exterior  hydrogen  stream 
expansion,  an  effect  of  the  gas-liquid  interaction  should  be  noted. 
Momentum  exchange  between  the  liquid  oxygen  spray  droplets  and  hydrogen- 
rich  gases  surrounding  them  accelerates  the  droplets  (Eq.  8a)  and 
decelerates  the  attendant  gases.  An  increasing  interior  cross-sectional 
area  is  needed  to  accommodate  the  decelerated  gas  flow;  the  surrounding 
hydrogen  stream  and  exterior  mixing  zone  must  be  displaced  radially 
outward  by  this  internal  expansion. 


In  the  model  formulation,  the  exterior  hydrogen  deceleration  was  assumed 
to  involve  increasing  mass  flowrate  of  hydrogen  by  having  the  velocity 
defect  penetrate  to  smaller  radius: 


rdef,2  '  rdef,l 


-c, 


Ax 


(13) 


def  ,2 
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a  value  of  from  free-jet  experimental  results  was  used.  The  increase 

across  the  Ax  increment  of  the  spray  zone  radius,  r  (whose  evaluation 

s  z 

is  discussed  in  a  following  section),  and  the  decrease  in  mass  flowrate 

of  hydrogen  still  having  injection  velocity,  w„  ,  were  then  used  to 

<2P 

calculate  the  effect  of  spray  zone  expansion  on  the  radius  to  which  the 
hydrogen  stream  velocity  defect  penetrates: 


def  ,2 


w, 


(rdef,l 


rszl> 


V 

HgPl 


+  r 


sz2 


1/2 


(15) 


The  maximum  hydrogen  stream  velocity  was  defined  as: 


t\) , 


U, 


V 


w, 


V2 


>  0 


pax 


U 


gb2 


w”  ,o  =  0 


V2 


(16) 


A  linear  approximation  of  the  Gaussian  velocity  distribution  in  the 
external  velocity  defect  flow  was  used  to  estimate  the  radius  of  the 
boundary  between  the  cylindrical ly  symmetric  propellant  flow  system  and 
the  surrounding  (stagnant)  chamber  gases: 
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Downstream  of  x  ,  the  furthest  distance  to  which  hydrogen  at  injection 
P 

velocity  penetrates,  continued  momentum  exchange  with  the 

spray  and  chamber  gases  must  result  in  deceleration  of  the  entire  gas/ 
spray  flow  system.  As  noted  earlier,  the  maximum  velocity  for  an 
unconfined  gaseous  jet  would  decay  inversely  with  increasing  distance. 
With  the  liquid  spray  core,  however,  such  a  rate  of  deceleration  would 
soon  result  in  gas  flow  velocities  lower  than  the  liquid  oxygen  spray 
droplet  velocities;  droplet  deceleration  would  then  transfer  momentum 
back  to  the  hydrogen-rich  gases  and  prevent  them  from  decelerating  so 


rapidly.  In  the  model  formulation,  therefore,  the  gas  velocity  at  the 


boundary  between  the  spray  zone  and  the  surrounding  flow  was  assumed  to 


follow: 


Values  of  the  exponent,  c^,  between  0.5  and  1.0  were  chosen  to  obtain 
deceleration  to  velocities  lower  than  a  calculated  turbulent  flame  speed 
upstream  of  the  position  where  the  abrupt  velocity  step  occurred  o.  the 
streak  photographs.  —  — 

A  conceptual,  schematic  cross-section  of  the  axial  propellant  velocities 
is  shown  in  Fig.  2.  The  solid  line  shows  the  concentric,  one-dimensional 
liquid  oxygen  and  gaseous  hydrogen  Hows  with  rectangular  velocity 
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profiles*  nt  the  injector  face.  A  short  distance,  x^,  downstream  some 

spray  has  been  formed,  resulting  in  internal  deceleration  of  hydrogen 

nnd  radial  displacement  of  the  whole  structure.  Somewhat  further  down= 

stream  at  x  .  the  unmixed  hydrogen  jet  has  just  been  engulfed  by  the 
P 

simultaneous  momentum  requirements  of  the  accelerating  spray  nnd  the 
external  gas  mixing  zone.  The  continued  deceleration  nnd  expansion  even 
further  downstream  is  shown  at  a  distance  from  the  injector. 


Rec  i  rcu  1  at  i oil 


Much  as  the  classical  bourded  circular  gns  jet  induces  an  annular  eddy 
of  gases  circulated  upstream  along  the  bounding  walls  and  ingested  into 
a  mixing  zone  around  the  jet,  chamber  gas  circulation  must  be  induced  by 
the  coaxial  jet  gaseous  hydrogen  injection.  Upstream  circulation  and 
entrainment  of  combustion  product  gases  in  this  manner  is  usually  referred 
to  as  recirculation. 

Just  what  gases  are  circulated  upstream,  their  composition,  and  tempera¬ 
ture  may  be  critically  important  to  the  steady-state  propellant 
combustion  nnd  precombustion  processes.  Recirculated  combustion  proJuct 
gas  would  provide,  ns  a  portion  of  its  mass,  additional  gaseous-state 
oxidant  to  the  coaxial  propellant  system  nnd  would  enhance,  through  its 
high- temperature  thermal  content,  heat  transfer  from  the  hydrogen  stream 
mixing  zone  to  the  liquid  oxygen  spray.  This,  in  turn,  would  enhance 
thi  spray  evaporation  rate;  conce  i  vn*»  ly ,  the  propellants’  attainment  of 
combustible  conditions  could  be  greatly  hastened. 
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Conversely,  circulation  of  uncombusted  gases  could  play  little  more 
than  a  passive  role  in  the  combustion  chamber  processes.  Such  a 
circulation  was  hypothesized  from  the  experimental  observations  (Ref.  1) 
and  is  reflected  in  the  schematic  representation  of  Fig.  1.  If  the  re¬ 
circulating  gases  in  the  annular  eddy  are  pure  hydrogen,  they  can  have 
no  appreciable  effect  on  the  propellant  precombustion  processes.  If, 
however,  nil  of  the  hydrogen  stream  were  drawn  into  the  spray  zone  up¬ 
stream  of  the  end  of  the  recirculation  eddy,  the  circulating  gases  would 
contain  some  oxygen;  a  moderately  higher  gaseous  oxygen  concentration 
throughout  the  gaseous  mixing  zone  and  the  spray  zone  (between  x  and 
the  downstream  end  of  the  eddy)  should  then  be  found. 

The  current  model  was  formulated  without  accounting  for  either  uncombusted 

gas  circulation  or  combustion  gas  recirculation.  The  chamber  volume 

outside  of  r  (Eq.  17)  was  simply  assumed  to  be  filled  with  stagnant 
"o~ 

gases  of  unspecified  nature  which  decelerated  a  "mixing  zone"  portion  of 
the  hydrogen  stream  without  actually  getting  mixed  with  the  hydrogen. 

These  assumptions  (synonymous  with  a  constant-mass,  constant-momentum 
eddy)  were  unavoidable  in  a  model  based  on  a  stepwise  progression  down 
the  chamber  length.  The  assumed  structure  is  very  nearly  realistic  if 
pure  hydrogen  gases  are  circulating. 

A  criterion  for  determining  whether  or  not  recirculation  forms  an 
important  part  of  the  flow  processes  associated  with  a  dected  turbulent 
gas  jet*  has  been  developed  by  Becker,  et  al.  (Ref.  10).  They  defined 


*A  ducted  jet,  in  general,  consists  of  a  gas  stream  flowing  through  a 
duct(ground  stream)  with  a  gas  jet  having  higher  or  lower  velocity 
injected  into  it,  usually  in  parallel  flow.  The  bounded  or  confined  jet 
corresponds  to  the  case  of  no  ground  stream  flow. 
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a  similarity  parameter,  Ct,  called  the  Crnyn-Curtet  number  and  determined 
experimentally  that  recirculation  occurred  in  n  constant-pressure, 
isothermal,  single-component  system  only  if  the  value  of  Ct  was  less  than 
3  4.  The  smaller  the  value  of  this  number,  the  stronger  was  the 
rec ireulati on. 

For  a  confined  gas  jet,  Ct  was  given  as: 


Ct  = 


"k((y*j  -.aH  v  j 


1  2 


(19) 


The  kinematic  velocity  (l'^)  in  a  cylindrical  duct  is  the  space  mean 

velocity,  which  is  related  through  mass  continuity  to  the  initial  jet 

o  <) 


velocity , 


l .  e .  ,  U  . 

I 

2  i »  oJ 

r  U  .** 

J  J 


r  .“  =  U.  r  Substituting  this,  together  with 
J  k  w 

into  Eq.  19  gives: 


Ct 


1 


1  -  1  2  ( 


r 

J 


-1  2 


(20) 


showing  that  recirculation  or  eddy  strength  for  a  confined  gas  jet  in  a 
cylindrical  duct  at  constant  pressure  is  dependent  only  upon  system 
geome  t  ry . 


The  two-phase  two-component  coaxial  injection  under  consideration  here 
differs  from  a  confined  gas  jet  mainly  in  that  momentum  is  transferred 
between  the  gas  stream  and  the  liquid  phase.  The  gas-phase  momentum 
available  for  inducing  circulation  is  thus  a  decaying  function  of 
distance  from  the  injector.  Becker's  result  expressed  in  Eq.  19  is  here 
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modified  by  replacing  I  with  (h  I  ),,.  Taking  the  ratio 

J  K 

(T  I  )  I  =  F  .  the  Craya-Curtet  number,  in  a  form  related  to 
v  g72  Ji  12’ 

Eq.  20  is: 


1/2 


Ct0  = 


rll02  _  rlj‘2 


v2 


12 


2  r 


-  1  2 


(21) 


Ihe  model  calculations  included  computation  of  this  value  of  the  Craya- 
Curtet  number  at  each  axial  position.  It  was  used  simply  as  an  indicator 
of  the  model's  validity:  If  flammable  conditions  were  reached  after 
Ct()  exceeded  3  4 ,  the  omission  of  circulation  from  the  model  was  assumed 
to  be  a  valid  simplification.  On  the  other  hand,  if  flammable  conditions 
were  calculated  to  occur  before  Ct,,  exceeded  3  4  ,  this  was  interpreted 
ns  a  warning  that  combustion  gas  recirculation  wns  probable  and  that 
the  model  calculations  might  be  grossly  in  error. 


Liquid  Oxygen  Spray  Zone 

The  structure  of  the  liquid  oxygen  spray  zone  surrounding  each 
cylindrical  liquid  oxygen  jet  is  not  known  with  any  certainty.  Attempts 
to  determine  the  spray  zone  diameter  ns  a  function  of  distance  from  the 
injector  (by  direct  measurements  on  the  high-speed  Fastax  motion  picture 
films)  gave  fairly  erratic  results;  both  jet-to-jet  diameter  variations 
and  f rame-to- f rame  variations  on  an  individual  jet  spray  system  contrib¬ 
uted  to  wide  measurement  scatter.  Further,  it  is  probable  that  sparse 
droplet  populations  in  the  outermost  periphery  of  the  spray  were  not 
visible  in  the  photographs.  Thus,  only  general  guidelines  for  the  spray 
zone  formulation  could  be  derived  from  the  experimental  work. 
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Certa  nly  the  gas  velocity  in  the  spray  zone  is  highest  at  (or  near) 
its  outer  periphery  and  is  lowest  immediately  next  to  the  liquid  oxygen 
jet.  This  is  shown  in  Fig.  11  as  a  continuously  decreasing  velocity,  from 
boundary  velocity  .it  r^  to  liquid  oxygen  injection  velocity  at  r^ 

The  radial  velocity  profile  was  hypothesized  in  the  model  formulation  to 
be  similar  to  that  across  a  boundary  layer.  Since  no  experimental  or 
theoretical  information  was  found  in  the  literature  for  a  comparable 
heterogeneous  turbulent  boundary  layer,  it  was  reasoned  that  the  boundary- 
layer-  thi  eke.. .  ve  loc  1  ty-gradient-reduc  ing  effects  of  spray  in  the 
gases  would  be  oi.  j  partially  offset  by  the  opposite  effects  of  having 
cylindrical  flow.  Consequently,  the  spray  zone  velocity  profile  was 
assumed  to  be  similar  to  that  for  a  laminar  boundary  layer  on  a  flat 
plate  ( Hef .  11): 


Conceivably,  a  sufficiently  realistic  description  of  the  spray  zone 
might  result  from  assumptions  that:  (l)  spray  zone  gases  are  homogeneous, 
(j)  all  sizes  of  pray  droplets  are  uniformly  mixed  with  the  homogeneous 
gas  flow,  and  (l)  liquid  oxygen  jet  atomization  and  spray  droplet  accel¬ 
eration  can  be  calculated  using  an  appropriate  mean  gas  velocity.  After 
an  initial  attempt  to  base  a  model  formulation  on  these  simplifications*, 


*The  preceding  assumptions  were  made  together  with  two  other  assumptions: 
that  the  spray  could  be  considered  as  monodisperse  and  that  the  spray 
zone  radius  could  be  prescribed  by  measuremei.  from  the  high-apeed  motion 
pictures.  The  calculated  results  showed  this  approach  to  be  so  oversim¬ 
plified  that  the  entire  line  of  attack  was  abandoned.  By  far,  the  poorest 
of  the  simplifications  was  the  use  of  a  monodisperse  spray;  probably  the 
next  poorest  was  basing  acceleration  and  evaporation  on  a  me  n  gas  velocity. 
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it  was  feared  that  too  many  details  of  the  spray  evaporation,  gaseous 
oxygen  distribution,  and  hydrogen  deceleration  processes  ver**  being 
averaged  out.  A  more  complex  model  of  the  spray  zone  was  thus  adopted. 


The  spray  zone  structure  in  the  final  formulation  was  hypothesized 
to  consist  of  thin,  cy 1 i ndr i cal ly  concentric,  annular  "shells"  of  spray, 
the  new  sprn>  being  shorn  from  the  liquid  oxygen  jet  surf nee  in  the 

Ax-increment  from  the  injector  was  considered  to  occupy  an  nnnulnr 
space  immediately  adjacent  to  and  surrounding  the  liquid  oxygen  jet 
(Fig.  3a).  To  make  room  for  this  new  spray,  the  spray  formed  in  the 
(k-1)^  Ax-increment  (upstream)  had  to  have  been  moved  out  radially 
into  a  larger-diame ter  annulus.  Successively  older  spray  elements  were 
presumed  to  be  flowing  through  annuli  lying  further  and  further  radially 
outward  from  the  liquid  jet  with  the  earliest-formed  spray  occupying  the 
outer  periphery  of  the  spray  zone. 


Consider  the  spray  element  being  formed  in  the  k 


th 


Ax-increment:  Sinc» 


it  is  new  sprny,  there  is  no  gas  flow  associated  with  it;  the  nearest 
available  gas  is  that  flowing  in  the  (k-l)*^  spray  element.  The  energy 
for  atomizing  and  the  momentum  for  accelerating  the  k^*  spray  group 
must  be  derived  from  influx  of  some  (or  all,  or  even  more  than  all)  of 
the  gas  previously  surrounding  the  (k-l)**^  sprny  droplets.  That  gas  must 
then  surround  the  k^*1  sprny  element  droplets  to  provide  sprny  momentum; 
in  so  doing,  the  composition,  temperature,  mid  velocity  of  the  gas  are 
changed  by  its  association  with  that  new  spray. 


At  the  same  time,  the  (k-l)^  spray  droplets  must  be  displaced  radinllv 

outward  by  the  deceleration  and  expansion  of  their  former  gases  thnt 

are  being  transferred  into  the  k^*1  sprny  group.  But  Eq.  22  states 

that  the  gas  velocity  is  higher  at  the  increased  spray  zone  rndius,  so 
t»H 

the  (k-l)  droplets  will  be  accelerated.  An  influx  of  gas  from  the  (k-2; 
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sprny  element  is  thus  required,  both  to  make  up  for  the  gas  transferred 
inwardly  and  to  provide  the  increased  liquid  oxygen  spray  momentum. 
Continuation  of  this  conceptual  process  leads  to  the  already  obvious 
conclusion  that  the  outermost  spray  groups  require  an  influx  of  hydrogen 
to  satisfy  their  gas  mass  and  momentum  requirements. 

Spray  /.one  structure  was  computed  by  means  of  a  complex  portion  of  the 
digital  computer  program.  The  programmed  equations  constitute  a  detailed 
system  for  ensuring  continuity  of  mass  and  energy  during  the  in tere lemeut 
gas  exchange  processes.  Since  the  accounting  system  is  not  of  direct 
interest  to  the  combust  ion  processes  per  se ,  only  a  description  of  the 
conceptual  bases  of  the  spray  /.one  solution  need  be  developed  here.  This 
is  accomplished  in  the  next  few  paragraphs  by  including  this  portion  in 
a  complete  summari/at  m  of  the  computer  model  calculation  scheme  for 
the  nonburning  region. 


Summary  of  Model  Calculation  Scheme 

As  previously  noted,  the  digital  computer  program  is  based  on  using 
knowledge  of  conditions  at  a  plane  x^  to  calculate  the  conditions  at  a 
plane  x,,,  a  distance  Ax  downstream.  A  simple  stepwise  progression  is 
used;  the  value  of  A  is  j  esumed  to  be  small  enough  that  mathematical 
instabilities  and  accumulation  of  appreciable  errors  are  avoided. 

The  program  first  calulates  liquid  oxygen  jet  atomization,  if  the  maximum 
droplet  diameter  is  not  calculated  to  be  exceeded.  The  acceleration  (or 
deceleration)  of  each  group  of  sprny  droplets  is  then  computed,  followed 
by  calculation  of  each  group's  evaporation  rate.  Next,  the  momentum 
changes  of  the  spray  droplets  and  new  oxygen  vapor  are  summed  to 
obtain  an  estimate  of  how  much  hydrogen  must  be  decelerated  and  mixed 
into  the  spray  zone.  That  information  is  used  to  approximate  the  new  sprny 
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zone  radius  and  boundary  velocity  which  are  then  used  in  a  single 
iteration  of  the  procedure  to  obtain  an  improved  estimate  of  spray  zone 
radius . 

The  total  spray  zone  area  is  apportioned  among  the  individual  spray 
elements  on  a  basis  of  propoi tionate  combined  gas  and  spray  mass  flow- 
rate.  The  downstream  radii  of  each  spray  element  annulus  are  calculated, 
and  a  new  mean  gas  velocity  associated  with  each  element  is  then  computed 
using  its  mean  downstream  radius.  The  resultant  distribution  of  spray 
element  areas  and  gas  velocities  is  assumed  to  be  inviolably  determined 
by  the  downstream  spruy  zone  rudius  and  boundary  velocity. 

The  gas  exchange  and  mixing  processes  shown  in  Fig.  3  are  then  calculated, 
beginning  with  the  last-formed  spray  group  and  proceeding  radially  outward 
to  successively  older  spray  groups.  The  gas  flowing  out  of  an  element 
is  presumed  to  result  from  mixing  gases  from  three  sources:  (l)  that 
which  flowed  into  the  element  from  upstream,  depleted  by  the  amount 
required  by  inly'ing  elements,  (2)  evaporation  of  spray  flowing  through 
the  element,  and  (3)  gases  required  from  outlying  elements.  Each  cal¬ 
culation  proceeds  ns  if  the  first  two  quantities  are  invariant;  the  value 
of  the  third  is  adjusted  by  iteration  until  the  gns  mixture  flowrate, 
tempe rature ,  and  density  satisfy  the  prescribed  exit  velocity  and  area 
requirements.  Equations  providing  for  mixing  among  the  three  streams 
while  preserving  continuity  of  hydrogen  and  oxygen  mass  flow  and  energy 
were  required.  The  physical  nnd  thermodynamic  properties  needed  for  solving 
the  spray  element  mixing  equations  nre  detailed  in  Appendix  A. 

The  initial  conditions  for  the  spray  /.one  are  preserved,  so  that  if  the 
calculation  results  are  unsatisfactory,  an  appropriate  adjustment  can 
be  made  nnd  the  calculations  reperformed.  Two  types  of  readjustment 
are  very  commonly  encountered:  (l)  a  spray  element  may  need  to  acquire 
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additional  gas  flow  when  there  is  none  available  (successively  smaller 
\  allies  of  spray  zone  radius  Bust  therefore  be  assumed  and  the  spray 
element  areas,  velocities,  and  gas  exchanges  be  recomputed  until  the 
requirements  are  satisfied  and  mass  ba lances);  and  (2)  excessive  gas  flow 
mav  be  calculated  to  exist  in  the  spray  zone,  usually  because  of 
deceleration  to  gas  velocities  lower  than  spray  eler  it  velocities. 

In  this  case,  the  spray  zone  radius  is  increased  to  accomodate  more  gas 
at  the  lower  velocities,  and  the  spray  element  areas,  velocities,  and 
gas  exchanges  are  recomputed. 

After  a  satisfactory  solution  for  spray  zone  structure  has  been  achieved, 
the  external  hydrogen  stream  expansion  and  deceleration  are  computed. 
Changes  in  spray  droplet  temperatures  are  calculated.  Finally,  an 
overall  momentum  balance  is  used  to  calculate  the  change  of  static 
pressure  across  the  A x- increment . 

The  calculated  conditions  at  the  downstream  end  of  the  A x- i ncrement 
are  examined  for  fulfillment  of  combustibility  criteria.  If  the  criteria 
are  met,  the  nonburning  model  is  replaced  by  a  combustion  region  model; 
otherwise,  the  calculated  results  are  used  as  initial  conditions  for 
extending  the  nonburning  region  over  another  A x- increment . 

CONDITION?  FOR  COMBUSTIBILITY 

In  Ref.  1,  it  was  stated  that  ”at  least  two  criteria  must  be  satisfied 
before  a  flame  can  penetrate  the  spray  zone  and  combustion  can  be 
sustained:  first,  enough  oxygen  must  have  evaporated  so  that  the 

gases  within  the  spray  .’one  are  flammable  (perhaps  5  to  7  volume  percent 
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oxygen),  and  second,  the  axial  gas  velocity  must  be  equal  to  or  less 
than  (he  turbulent  flame  speed,  or  else  nn  established  tombusiio 
would  be  swept  downstream."  Further  discussion  developed  the  probability 
that  any  flame  resulting  from  ignition  in  the  interior  of  the  spray  zone 
is  not  a  steady- state  possibility  so  long  as  any  of  the  nonburning  spray 
zone  gas  is  flowing  at  velocities  in  excess  of  the  turbulent  flume  speed. 
This  led  to  the  suggestion  that  a  third  condition  for  combustibility 
might  be  the  existence  of  a  strong  stabilized  flame  source,  such  as  a 
continuous  eddy  of  recirculating  combustion  gases  surrounding  the  spray 
/one . 

These  three  conditions  were  used  in  the  model  formulation  for  determining 
combustibility  of  the  spray  zone  gases.  The  tlird  condition  was  satisfied 
by  requiring  that  the  ges  velocity  in  the  outermost  spray  element  which 
still  had  unevapovated  liquid  oxygen  sprnv  in  it  not  be  greater  than  the 
calculated  turbulent  flame  velocity.  If  the  gas  velocities  throughout 
the  sproy  zone  were  below  the  turbulent  flame  speed,  combustibility  was 
considered  to  be  satisfied  when  the  oxygen  concentration  in  any  spray  zone 
element  satisfied  the  flammability  limit  criterion. 

The  values  for  the  upper  flammability  limit  for  gaseous  oxygen  and 
hydrogen  mixtures  used  in  the  model  calculations  rnnged  from  5  to  8 
\olume  (or  mole)  percent  oxygen.  The  turbulent  flame  velocity  was 
assumed  to  be  dependent  only  upon  pressure: 

Ut  =  a  P1  2  (23) 

A  detailed  discussion  of  the  selection  of  these  relationships  is  given 
in  Appendix  C. 
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At  that  position  in  the  combustion  chamber  where  the  spray  zone  com¬ 
bustibility  requirements  are  satisfied,  the  nonburning  injection  region 
model  is  replaced  by  a  combustion  region  model.  Some  general  features 
of  this  region's  structure  are  discussed  briefly  before  the  details  of 
the  model  formulation  are  developed. 

Initial  conditions  for  the  combustion  region  are  established,  and  these 
are  used  to  calculate  conditions  one  Ax-increment  further  down  the 
chamber.  The  model  calculations  thus  continue  stepwise  along  the  length 
of  the  chamber  in  the  manner  previously  discussed. 


Flame  Spread i Tig 

Since  the  axial  gas  velocities  within  the  spray  zone  must  be  lower  than 
the  turbulent  flame  speed,  flame  should  propagate  from  the  point  of 
ignition  across  the  entire  spray  zone  within  an  axial  "distance  that  is 
of  the  order  of  the  spray  zone  thickness.  With  injection  element  spacings 
on  the  order  of  1/2  inch,  spray  zone  thicknesses  of  the  order  of  l/8  to 
1 A  inch  appear  reasonable.  The  chamber  length  associated  with  flame 
spreading  through  the  spray  zone  is  thus  only  a  very  small  fraction  of 
the  chamber  length;  simplification  to  the  assumption  of  a  plane  flame 
front  across  the  spray  zone  at  the  position  of  first  combustibility  thus 
appears  to  be  justified.* 


*Some  model  calculations  were  performed  in  which  ignition  was  provided 
for  only  in  the  gases  directly  against  the  liquid  orygen  jet  and,  even 
though  the  oxygen  concentrations  in  the  outer  elements  were  too  low  to 
satisfy  flammability,  expansion  of  hot  gases  out  of  the  core  of  the  spray 
zone  increased  the  spray  gasification  rates  such  that  the  combustible 
region  spread  radially  outward  as  rapidly  as  the  more  likely  inward 
propagati on. 
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Spray  Zone  Combu s t i on 


Superposition  of  combustion  on  the  annular-shell  structure  of  the 
nonburning  spray  /.one  appeared  to  impose  insurmountable  difficulties. 

The  heterogeneous  combustion  probably  generates  intense  local  turbulence 
(hef.  12)  tlmt,  would  invalidate  such  an  analysis  in  any  case.  The 
contents  of  the  nonburning  spray  zone  elements,  therefore,  were  considered 
to  be  thoroughly  mixed  upon  passage  through  the  standing  plane  flame 
front.  The  mixed  gases  were  presumed  to  burn,  resulting  in  a  one- 
dirmnsional  flow  of  combustion  gases  with  all  of  the  sprays  uniformly 
mixed  in  it. 


The  combustion  process  was  accounted  for  in  an  indirect  way  similar  to 
that  used  for  liquid  oxygen/ilP-1  propellants  earlier  (Ref.  5).  At  each 
local  axial  plane  across  the  spray  combustion  zone,  the  flowrates  of 
gaseous  hydrogen  and  gaseous  oxygen  that  had  been  mixed~to  form  the 
combustion  gas  stream  were  considered  to  have  reacted  to  thermodynamic 
equilibrium,  stagnation  conditions,  and  then  to  have  been  accelerated 
isentropically  to  the  local  flow  conditions.  This  permitted  considerable 
simplification  of  the  combustion  model  in  that  all  of  the  desired 
combustion  gas  properties  could  be  tabulated  (or  curve-fitted  with 
polynomials)  as  simple  functions  of  local  gaseous  propellant  mixture 
ratio.  The  results  of  complete  thermochemical  performance  calculation 
programs  are  thug  used  without  introducing  the  complexity  of  those 
calculations  into  the  combustion  model.  In  Appendix  D,  values  of  the 
combustion  gas  properties  so  used  are  given.  Other  physical,  thermodynamic 
and  transport  properties  required  for  the  solution  of  the  combustion 
model  are  also  given  in  Appendix  D. 
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Transverse  I nhomogene  t  1 1  ph 

The  combustion  region  model  tl  on  must  account  for  three  concentric  flows: 
(1)  any  residual  liquid  oxygen  jet  that  penetrates  as  far  downst ream  as 
the  combustion  region,  (2)  the  spray  combustion  /one,  and  (3)  any  residual 
hydrogen  stream  or  other  gas  flow  which  surrounds  the  spray  combustion 
zone  but  does  not  itself  contain  liquid  oxygen  spray.  The  model  formul¬ 
ation  consists  of  expressions  for  the  interaction  and  mass  transfer 
between  the  spray  combustion  /.one  and  the  liquid  oxygen  jet  and  between 
the  spray  combustion  zone  and  the  surrounding  gas  flow,  ns  well  as  for 
the  internal  behavior  of  t lie  latter  two  streams. 

In  general,  two  kinds  of  inhomogeneity  might  be  expected  to  occur  and 
persist  for  some  distance  axially  into  the  combustion  region.  The  first 
is  implied  by  recognition  of  a  gas  stream  that  does  not  contain  spray. 
Initially,  this  stream  lay  have  been  only  slightly  more  fuel  rich  than 
the  gases  in  the  spray  combustion  zone.  Further  downstream,  however,  its 
fixed  composit  ion  will  become  more  and  more  fuel  rich  with  respect  to 
spray  combustion  /.one  gases  because  of  continued  oxygen  spray  gasification. 

The  second  kind  of  transverse  inhomogeneity  must  develop  within  the 
spray  combustion  /one  as  a  result  of  accelerated  atomization  and  gasi- 
f  nation  of  the  residual  liquid  jet.  The  atomization  equations,  (Eq.  la 
and  2a),  show  that  a  moderate  velocity  increase  resulting  from  burning 
the  gases  suiioundtng  the  liquid  oxygen  jet  will  both  increase  the  jet 
atomization  rate  and  decrease  the  mean  droplet  size  being  produced. 

Further,  the  gases  surrounding  the  spray  droplets  are  at  elevated 
temperature  so  that  increased  heat  transfer  to  the  drops  hastens  their 
gas i f i ( at i on .  A  residual  liquid  oxygen  jet  containing  n  large  portion 
of  the  injection  flowrate  may  thus  be  atomized  and  gasified  fairly 
rapidly  in  only  a  short  segment  of  chamber  length  to  form  a  highly  oxygen- 
ri<h  (ore  in  the  spray  combustion  zone. 
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Roth  of  those  kinds  of  inhomogene 1 ty  in  the  propellant  combustion  field 
will  be  degraded,  with  increasing  distance  from  the  injector,  by 
turbulent  mixing  processes.  Whether  a  completely  mixed,  fully  burned 
equilibrium  combustion  chamber  exhaust  gas  (which  is  desirable  for 
maximizing  performance  efficiency)  is  achieved  must  ultimately  be  deter¬ 
mined  by  this  type  of  gas-phase  mixing. 

In  formulating  a  model  of  the  combustion  region,  an  accounting  for  both 
kinds  of  inhomogeneitv  appeared  to  be  inordinately  complex,  if  not 
intractable.  The  further  simplification  was  introduced  that  the  oxygen 
spray  and  vapor  resulting  from  liquid  oxygen  jet  atomization  i9  continually 
mixed  uniformly  with  the  other  spray  and  combustion  gases  in  the  spray 
combustion  zone.  Only  turbulent  mixing  of  the  surrounding  fuel-rich 
gas  stream  into  that  oxygen-rich  zone  then  was  considered.  As  a  result 
of  this  simplification,  it  was  necessary  to  provide  for  accentuating  the 
apparent  effect  of  the  latter  process  by  diverting  some  of  the  nonburning 
spray  zone  gases  around  the  flame  front  und  into  the  surrounding  gas 
stream.  While  some  such  process  m  ;y  actually  occur,  its  accentuation 
here  should  be  considered  as  an  artifice  to  compensate  for  the  simplifica¬ 
tion  concerning  the  absence  of  an  oxygen-rich  core  in  the  spray  combustion 
zone . 


Initial  Conditions 


The  three  concentric  streams  referred  to  previously  were  initialized  at 
the  end  of  the  nonburning  region.  The  liquid  oxygen  jet  flows  from  one 
region  into  the  next  without  change.  Similarly,  the  liquid  oxygen  spray 
lements  are  unchanged  with  respect  to  sizes  and  numbers  of  droplets, 


Hit  64-20 


31 


A  DIVISION  OF  NORTH  AMERICAN  AVIATION 


I  N  C 


but  they  are  no  longer  presumed  to  be  separated  into  annular  strata. 
Before  initializing  the  spray  combustion  zone  gas  flow,  a  portion  of 
the  gas  of  each  nonburning  spray  element  was  diverted  radially  outward 
into  the  surrounding  gas  stream.  The  proportion  of  gas  flow  to  be 
diverted  out  of  the  nth  spray  element  was  arbitrarily  taken  to  be: 
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d  i  v  ,  n 
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k  +  1  -  n 
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where  k  represents  the  number  of  spray  elements  that  still  contain 
unevaporated  spray,  with  the  k^  element  lying  next  to  the  liquid 
oxygen  jet.  The  flowrate  of  gases  into  the  spray  combustion  zone  is 
then : 


w 

cgl 


(1  -  F 


di v ,n) 


and  the  amount  of  oxygen  in  that  gas  is: 


k 

*0,,cgl  ^  X0ognl  Wgnl  ^ 1  '  Fdiv,n^  Mgnl 

*•  n  =1 

By  simple  mass  cortinuity  relationships,  the  flowrate  and  oxygen  content 
of  the  surrounding  gas  stream  are  then  determined. 
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The  actual  combustion  of  the  spray  combustion  zone  gases  is  considered 
to  take  place  across  the  first  combustion  region  Ax- increment  so  that 
the  initial  conditions  are  still  nonburning.  The  static  pressure, 
therefore,  is  most  likely  the  same  as  that  last  calculated  for  the 
nonburning  region.  Continuity  of  momentum  then  requires  that  the  spray 
combustion  zone's  initial  gas  velocity  be: 
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cgl 


n-1 


(1 


F.  ) 
di v ,n' 


(-7 ) 


The  cross-sectional  area  required  for  this  gas  flow  is  then: 


(28) 


where  PCg^  is  obtained  from  a  perfect  gas  law  relationship 


144  P1 

R  T 


cgl 


(29) 


requiring  that  additional  mixing  equations,  paralleling  Eq.  26,  be  employed 
to  obtain  T  j  through  thermal  continuity. 

Another  set  of  equations,  similar  to  those  used  for  the  spr  y  combustion 
zone  gases,  was  written  for  simultaneously  obtaining  initial  values  of 
the  surrounding  gas  stream  flowrate,  composition,  temperature,  density, 
velocity,  and  flow  area. 
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If  the  sum  of  the  liquid  oxygen  jet,  spray  combustion  zone,  and  gas 
stream  flow  areas  thus  calculated  exceeded  the  total  cross-sectional 
area  available,  the  two  gaseous  flow  areas  ./ere  reduced  proportionately 
to  make  that  sum  equal  to  the  chamber  area  per  injection  element.  An 
iterative  procedure  was  then  used  to  lower  the  static  pressure  and  in¬ 
crease  the  gas  velocities  until  overall  mass  and  momentum  continuity 
relationships  were  satisfied. 


Liquid  Oxygen  Jet  Atomization 

The  equations  presented  previously  for  liquid  oxygen  jet  atomization 
late  and  local  mean  droplet  sizes  produced  in  the  nonburning  region 
were  used  .is  well  in  the  combustion  region.  The  one-dimensional 

I  t 

combustion  gas  velocity  aid  density  replnced  and  in  Eq.  la 

and  2a.  The  previously  imposed  maximum  droplet  diameter  was  also 
honored  in  the  combustion  region,  so  that  atomization  could  continue 
(or  be  resumed)  within  a  particular  Ax- 1  ncrement  only  if  droplets  smaller 

than  that  maximum  were  to  be  formed. 

In  the  Fortran  formulation  of  the  computer  program,  storage  locations 
v.ere  reserved  for  a  maximum  of  100  iniividual  spray  groups.  This  was 
believed  to  be  a  sufficient  number  for  most  practical  systems  since  it 
would  account,  for  example,  for  liquid  jet  penetrution  to  a  minimum  of 
5  inches  from  the  injector  if  Ax  were  0.05  inch,  or  to  10  inches  if 
Ax  were  0.10  inch.  If  the  conditions  or  controls  had  been  chosen  such 
that  some  anatomized  liquid  oxygen  jet  were  computed  to  survive  past 
the  formation  of  100  spray  groups,  all  of  that  remaining  jet  was 
dumped  into  the  100th  spray  group;  while  this  is  certainly  unrealistic, 
it  permitted  the  machine  calculations  to  continue  without  using  values 
foi  variables  taken  from  storage  corresponding  to  entirely  different 
van  able s . 
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Individual  Droplet  Behavior 

The  dynamic  behavior  oi"  all  the  droplets  in  each  size  group,  during  their 
passage  through  a  given  Ax- i ncrement ,  was  calculated  using  equations 
related  to  Eq.  6  through  12.  Droplet  ballistics  were  determined  directly 
from  Eq.  6  through  9  using  combustion  gas  properties  wherever  the 
subscript  "gnl"  occurs  there.  Similarly,  changes  in  droplet  bulk  temper¬ 
ature  were  evaluated  by  means  of  Eq .  10  and  11.  Droplet  mass  evaporation 
rates  were  calculated  by  an  equation  related  to  Eq.  12  but  derived  for 
high  evaporation  rates  into  a  gas  stream  that  does  not  contuin  any  of 
the  evaporating  species  (Ref.  6): 
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dn2  dnl 
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Ddnl  NuMnl  P1  AX 
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(30) 


Noting  from  Eq.  29  that  the  terms  p  g  ( 1 44  M^)  in  Eq.  12  and  (RT^) 
in  Eq.  30  are  equivalent  for  a  perfect  gas  reveals  that  these  two 
evaporation  equations  differ  only  in  their  terminal  bracketed  terms. 


If  the  local  static  pressure  exceeds  the  critical  pressure  of  oxygen, 
it  has  been  predicted  that  nearly  all  or  the  heat  transferred  to  the 
liquid  oxygen  spray  droplets  will  be  absorbed  as  sensible  heat  of  the 
liquid  rather  than  being  returned  to  the  gas  stream  as  latent  heat  of 
vaporization  (Ref.  0  and  13).  Droplets  entering  (or  formed  within)  the 
spray  combustion  zone  thus  may  be  gasified  by  transient  heating 
through  the  critical  temperature  rather  than  by  quas i steady-state  surface 
evaporation.  The  process  was  included  in  the  model  formulation  by 
abruptly  adding  to  the  combustion  gas  stream  the  mass  flowrate  of  any 
spiny  group  calculated  to  have  droplet  temperatures  exceeding  the  critic.  1 
temperature  of  oxygen. 
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Turbulent  Gas  Mixing 

In  addition  to  the  liquid  oxygen  spray  gasification,  the  combustion  gf  s 
flowrate  may  be  increased  by  admixture  of  gases  from  the  surrounding  gas 
stream.  A  rather  gross  approximation  of  the  turbulent  mixing  between 
these  streams  was  employed  in  the  model  formulation:  a  total  turbilent 
mixing  length,  L  ,  was  selected  by  examination  of  the  velocity  vs  chamber 
length  plots  of  Ref.  1,  and  the  surrounding  gas  stream  flowrate  was 
diminished  linearly  with  distance  over  that  length.  The  combustion  gas 
flowrate  was  thus: 
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and  the  surrounding  gas  stream  flowrate  was; 

(32 

This  type  of  mixing  was  presumed  not  to  occur  unless  the  propellant  flows 
filled  the  available  chamber  cross  section. 


gs2 


=  w 


(w  ) . 

'  gs  1  '  1 


Ax 


gsl 


36 


RH  6>i -2‘) 


A  DIVISION  OF  NORTH  AMERICAN  AVIATION 


I  N  C 


ROCKUTDVmL 


Spray  Zone  Combustion 

It  was  noted  earlier  that  propellant  combustion  was  accounted  foi  in 
the  model  formulation  by  considering  the  local  propellant  flowrates  to 
have  reucted  locally'  to  thermodynamic  equilibrium  stagnation  conditions, 
dependent  (for  a  given  range  of  chamber  pressures)  only  upon  the  local 
gaseous  propellant  mixture  ratio.  This  requires  knowledge  of  the  equi¬ 
valent  flowrates  of  oxygen  and  hydrogen  that  resulted  in  the  combustion 
gas  stream: 
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The  stagnation  properties  of  the  combustion  gas  were  then  found  ns 
detailed  in  Appendix  D. 
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Static  conditions  for  the  flowing  combustion  gas  obtained  from  lsentropic 
flow  relationships: 


(37) 


(38) 


38 


HR  64-29 


A  OIVISION  OF  NORTH  AMERICAN  AVIATION 


I  N  C 


ROCKETDVWB  • 


Combustion  of  the  Surrounding  Gas  Stream 

The  initial  oxygen  concentration  in  the  gas  stream  surrounding  the  spray 
combustion  zone  was  compared  to  the  upper  flammability  limit  for  deter¬ 
mining  whether  this  stream  would  be  burned  or  remain  unburned.  If  the 
stream  were  incombustible,  the  ga.»  stream  stagnation  temperature  and 
stagnation  pressure  were  assumed  to  remain  constant  and  equal  to  their 
initial  values  at  all  downstream  positions.  This  is  equivalent  to  ignor¬ 
ing  heat  transfer  into  the  gas  stream  as  n  flows  axially  along  the 
chamber . 

For  the  burned  gas  stream,  the  stagnation  pressure  was  taken  to  be  the 
initial  value,  but  the  stagnation  temperature,  stagnation  velocity  of 
sound,  molecular  weight,  and  ratio  of  specific  heats  were  determined  at 
tne  beginning  of  the  combustion  region  in  the  same  way  the  spray  zone 
combustion  gas  properties  were  obtained.  Since  no  source  was  provided 
for  mass  transfer  into  this  stream,  the  values  for  those  variables 
determined  upon  first  combuution  were  assumed  to  remain  invariant  through¬ 
out  the  succeeding  flow. 


Area  Distribution,  Gas  Velocities,  and  Chamber  Pressure 


Until  the  total  available  cross  section  of  the  chamber  has  become  filled 
with  flowing  propellants  and  combustion  gases,  gas  generation  and  heat 
release  can  be  absorbed  by  transverse  expansion  without  large  gas 
velocity  or  pressure  changes  occurring.  Although  equations  were  included 
in  the  model  formu.'ation  for  computing  the  progress  of  such  combustion, 
only  a  short  segment  of  a  conventional  rocket  combustor's  chamber  length 
is  needed  to  achieve  full  flow,  so  the  details  are  omitted. 
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Vhen  the  chamber  has  become  filled  with  axial  propellant  and  combustion 
ffns  flow,  continued  gasification  of  liquid  oxygen  and  combustion  reactions 
between  unlike  gas  streams  must  necessarily  refult  in  acceleration  of 
thf  flows  and,  simultaneously,  in  lowering  of  the  static  pressure.  The 
cross-sectional  areas  and  velocities  of  the  two  gaseous  streams  and  the 
pressure  are  all  unknowns  and  are  all  interrelated. 


The  c chous ti on  region  Model  solution  was  based  on  approximating  the 
distribution  of  the  chamber  cross-sectional  area  among  the  liquid 
oxyK"n  J c *  the  spray  combustion  zone,  and  the  surrounding  gas  stream 
by  assuming  initially  that  the  pressure  and  flow  Mach  numbers  across  the 
Ax-increment  are  constant.  Estimates  of  downstream  gas  densities  and 
velocities  were  thus  obtained  and,  from  mass  continuity,  the  required 
cross-sectional  flow  areas  calculated.  In  general,  the  sum  of  those 
calculated  areas  should  exceed  the  total  chamber  area,  so  the  two 
gaseous  flow  areas  were  reduced  proportionately: 
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These  approximations  were  used  as  starting  values  which  were  continuously 
modified  during  a  subsequent  iterative  calculation  of  the  two  gas 
stream  densities,  velocities,  and  chamber  pressure  at  plnne  2. 
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The  iterative  procedure  began  with  an  estimate  that: 
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which  results  in  a  total  propellant  momentum  of: 
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Therefore 
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The  pa-  stream  flow  has  been  assumed  to  be  lsentropic.  An  inversion 
of  the  lsentropic  pressure  relationship  gives: 
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and  finally: 
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All  these  newly  calculated  values  of  combustion  and  surrounding  gas 
stream  properties  were  then  used  to  re-enter  th?  system  of  equations  at 
Eq.  44  to  obtain  modified  values  of  the  propellant  momentum  sum,  chamber 
pressure,  etc.  Experience  showed  that  the  solution  usually  converged  to 
constant  values  after  four  or  five  iterations;  the  formulated  computer 
program  employed  seven  iterations. 

Values  for  all  parameters  of  interest,  having  been  calculated  across 
an  nurement  of  chamber  length,  they  were  written  out  and  then  used  as 
initial  values  for  calculations  across  the  next  increment.  This  process 
was  continued  until  either  the  throat  position  in  the  combustion  chumber 
exhaust  nozzle  was  reached  or  the  combustion  gas  velocity  exceeded  the 
local  sound  velocity.  While  these  two  criteria  for  terminating  the  model 
calculations  should  actually  be  synonymous,  there  are  valid  reasons  that 
they  may  appear  to  be  different.  First,  the  model  computations  start 


1U1  64-29 


43 


A  DIVISION  OF  NORTH  AMERICAN  AVIATION 


I  N  C 


ROCKETDVRI  *1  • 


from  the  injector  face  with  prescribed  values  of  chamber  pressure,  pro¬ 
pellant  flowrates,  and  chamber  and  nozzle  configurations.  The  prescribed 
values  may  be  incompatible,  e.g.,  assumption  of  too  high  an  injection 
end  chamber  pressure  (or  too  low  a  flovrate  or  too  large  a  throat)  would 
result  in  calculation  of  subsonic  gas  velocity  at  the  geometric  throat. 
Conversely,  assumption  of  too  low  a  chumber  pressure,  etc.,  would  give 
sonic  gas  velocity  upstream  of  the  geometric  throat.  It  may  thus  be 
necessary  to  modify  the  prescribed  initial  values  for  the  two  criteria 
to  come  acceptably  close  to  occurring  at  the  same  axial  location. 


Combustion  Chamber  and  Nozzle  Configuration 

The  cross-sectional  area  per  injection  element  of  the  combustion  chamber 
and  nozzle  convergence  must  be  specified  as  a  function  of  distance  from 
the  injector.  For  the  transparent  two-dimensional  combustion  chamber 
employed  in  the  experimental  investigation  (Ref.  l),  these  areas  are 
somewhat  in  doubt  because  of  pyrex  liner  melting.  The  initial  pretest 
geometry  was  modified  ns  detailed  in  Appendix  E  to  provide  input  values 
for  the  combustion  model. 


kh 


HR  V i-20 


TDVIXE 


A  DIVISION  or  NORTH  A  H  t  R  I  C  A  N  AVIATION 


I  N  C 


CALCULATED  RESULTS  AND  DISCUSSION 


BEST-FIT  CASES 

The  selected  model  calculation  results,  which  have  been  found  to  most 
nearly  duplicate  the  experimental  information,  are  presented  and  discussed. 
The  reasons  for  the  particular  valuation  of  some  parameters  are  indicated. 


Initial  Model  Conditions 


The  majority  of  the  combustion  model  cases  calculated  have  been  with  ini¬ 
tial  conditions  duplicating  the  experimental  injection-end  conditions  of 
one  of  the  two  experiments  that  yielded  essentially  complete  experimental 
information,  i.e.,  tests  417  and  431.  The  initial  model  conditions  for 
these  two  cases  are  sunmarized  in  Table  1. 

Injection-end  chamber  pressure  was  not  measured  during  either  of  these 
tests;  the  values  found  to  result  in  model-calculated  pressures  that  match 
the  measured  pressures  further  downstream  are  given  in  the  table.  Fur¬ 
ther,  the  pressure  listed  for  test  417  is  100  psia  lower  than  that  re¬ 
ported  in  Ref.  1  because  the  model  calculations  corroborated  an  apparent 
discrepancy  between  measured  parameters  indicated  there. 


Values  Selected  for  Arbitrary  Constants 

The  chamber  length  increment  was  taken  as  A  x  =  0.03  inch  on  the  basis 
that  100  increments  (and,  therefore,  100  spray  groups)  would  permit  atomi¬ 
zation  to  proceed  over  a  minimum  of  3  inches  of  chamber  length. 
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The  atomization  constants  were  taken  as  B^  =  10,  =  0.15,  =  5*0  and 

c0  =  10.0  (tq.  la,  2a,  and  l).  This  value  of  1)  is  two  orders  of  magni- 

to  A 

tude  larger  than  that  recommended  by  Mayer  (lief.  2);  it  was  necessary  to 
avoid  forming  such  extremely  small  spray  droplets  ‘hat  their  acceleration 
would  result  in  very  abrupt  deceleration  of  the  hydrogen  stream.  Further, 
it  was  also  necessary'  to  avoid  such  rapid  evaporation  of  oxygen  that  com¬ 
bustibility  requirements  would  be  satisfied  very  near  the  injector.  Mayer 
indicated  that  in  Fq.  1  the  importance  of  liquid  viscosity  may  be  over¬ 
emphasized,  i.e.,  that  B^  may  depend  upon  ^  \  Because  the  liquid 

waxes  used  in  He f .  3  were  20  to  70  times  more  viscous  than  liquid  oxygen, 
this  effect  may  partially  account  for  the  very’  large  lalue  of  I)  required 
here.  As  noted  earlier,  some  additional  difference  might  be  expected 
because  Mayer's  evaluation  of  was  based  on  a  mean  droplet  diameter  for 
the  spray,  and  it  is  being  used  here  to  calculate  a  local  individual 

droplet  diameter.  The  value  of  C,  is  toward  the  high  side  of  the  range 

A 

recommended  in  Itef.  4. 


The  rate  of  decrease  of  the  apparent  radius  of  a  hydrogen  stream  core 
still  having  injection  velocity  was  controlled  by  taking  c^  =  0.11852 
in  Ixj.  13.  The  substquent  deceleration  of  the  boundary  between  the 
spray  zone  and  the  surrounding  flow  was  fixed  by  the  exponent  in  Eq .  18; 
a  value  of  c^  =  0.63  was  used.  As  noted  earlier,  this  exponent  would  be 
unity  for  a  free  jet;  a  value  lower  than  unity  seems  likely  to  result 
from  momentum  exchange  to  the  decelerating  gases  from  concurrently 
decelerating  liquid  spray. 


The  maximum  spray  droplet  size  was  selected  as  230  microns,  as  detailed 

in  Appendix  B.  Combustibility  conditions  were  taken  as  8.0  mole  percent 

1  2 

oxygen  and  =  12.0  P  (This  latter  equation  gives  flame  speeds 

approximately  40  percent  higher  than  the  turbulent  curve  in  Fig.  C-l . ) 
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The  fruitions  of  nonburning  spray  zone  gases  diverted  around  the  flame- 
front  and  into  the  gas  stream  surrounding  the  spray  combustion  zone  *Eq. 
24)  was  determined  by  c_,  -  0.35;  the  particular  value  was  selected  to 
limit  the  amplitude  of  the  abrupt  combustion  gas  velocity  jump  when  con¬ 
ditions  comparable  to  test  431  were  calculated.  The  total  turbulent 
mixing  length,  9  inches,  wus  also  (selected  primarily  from  the  experi¬ 

mental  \elocity  profile  of  that  test. 


Nonburning  Region 

Structure .  The  propellant  flow  structure  calculated  for  comparison  with 
test  417  ip  shown  in  Fig.  4.  (Note  that  the  radial  scale  ie  enlarged  by 
2-1  2  times  compared  to  the  axial  scale.)  As  may  be  inferred  from  the 
number  of  annulur  gus-spray  mixture  groups,  primury  utomization  of  the 
liquid  oxygen  jet  was  calculated  to  have  proceeded  over  approximately  the 
first  3/4  inch  from  the  injector  and  then  to  have  ceased.  The  first 
formed  spray  groups  were  calculated  to  have  been  accelerated  rapidly  and 
thus  to  have  decelerated  large  fractions  of  the  hydrogen  stream,  with  a 
resultant  radial  expansion  of  the  decelerated  gases.  In  fact,  the  en¬ 
tire  hydrogen  stream  flow  was  very  rapidly  engulfed  by,  or  drawn  into, 
the  expanding  spray  zone. 

It  is  interesting  to  compare  this  calculated  structure  with  nonburning 
experimental  observations  of  coaxial  jet  injection.  Figure  5  shows  two 
microflash  photographs  of  n  cylindrical  liquid  nitrogen  stream  injected 
into  umbier.t  uir  (lief.  14).  Only  the  3,'16-inch  diameter  l’quid  nitrogen 
jet  was  flowing  in  Fig.  5a »  while  an  annular  gaseous  helium  flow  (from 
a  325-psig  source  tank)  surrounded  the  liquid  nitrogen  jet  in  Fig.  5h.  The 
liner  atomization  imposed  by  the  helium  flow  is  apparent  and  the  corres¬ 
pondence  between  the  observed  outer  spray  boundary  there  and  the  calculated 
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structure  in  Fig.  4  is  rtriking.  It  was  noted  (Ref.  14)  that  an  unatomized 
liquid  core  appeared  to  persist  for  some  inches  further  downstream  with  the 
helium  flowing  than  without  it.  In  another  experiment,  two  parallel  cyl¬ 
indrical  jets  of  water,  spaced  0.53  inch  upart  and  surrounded  by  coaxial 1\ 
injected  gaseous  nitrogen,  were  observed  to  maintain  their  identities  as 
separate  liquid  jets  for  about  18  to  20  inches  before  jet  instability 
eventually  led  to  rather  coarse  atomization  (Fig.  b).  These  experimental 
observations  lend  su, oort  to  the  calculated  cessation  of  surface  shear 
atomization  a  short  distance  downstream  of  the  injector. 


Liquid  Oxygen-det  Atomization.  The  calculated  mean  sizes  of  spruy  drop¬ 
lets  produced  locally  by  the  surface-shear,  wave-crest-shedd  ng  mechanism 
(Fig.  7a)  remain  constant  as  long  us  some  hydrogen  stream  persists  at  in¬ 
jection  velocity;  it  becomes  larger  us  the  spray  zone  boundary  velocity 
decays  , 


The  distributions  of  spray  muss  produced  among  all  the  droplet  sizes  are 
shown  in  Fig.  7b.  The  initial  sizes  und  mass  atomization  rutes  were  used 
there  without  accounting  for  reductions  due  to  vaporization.  Figure  7b 
18  P1  otted  on  a  normalized  or  reduced  diameter  basis  to  facilitate  com¬ 
parisons  with  standard  distribution  functions.  The  three  meun  diameters 
noted  in  Fig.  7b  were  derived  from  the  calculated  droplet  masses  using 
the  relationship  (Ref.  15), 
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Because  the  volume- number  mean  diameter,  D^g,  6q®  been  found  to  compare 
well  with  experimentally  measured  mean  diameters  in  un  operating  rocket 
environment  (Ref.  16),  that  mean  diameter  was  used  for  normalizing  the 
calculated  droplet  diameters.  Ingebo  has  found  that  the  distributions  of 
spray  mass  from  impinging  doublet  (Ref.  17)  and  showerhead  injection  ele¬ 
ments  (Ref.  16)  were  well  described  by  the  Nukiyama-Tanasawa  distribution 
equat ion 


dv 

dD 


h6  6  D5 

r  (6  6) 


The  cumulative  volume  fraction  of  spray  smaller  than  any  arbitrary  dia¬ 
meter  is  showi  in  Ref.  1*5  to  be  given,  for  this  distribution,  by  the  ratio 
of  the  incomplete  and  complete  gamma  functions: 


v 


r  (6/6) 
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Evaluated  at  6  --  1,  this  expression  gave  the  indicated  Nukiyama-Tanasawa 
distribution  curve  in  Fig.  7b. 


It  is  seen  that,  because  the  relutive  velocit;  between  the  liquid  oxygen 
jet  and  the  hydrogen  stream  is  initially  very  high,  the  atomization  proc¬ 
ess  produces  a  spray  in  winch  both  the  numbers  of  droplets  and  the  spray 
mass  of  tiie  smaller  droplet  sizes  predominate.  While  most  rocket  injec¬ 
tors  produce  sprays  with  between  75  und  150  microns,  the  surface-wave 
mechanism  produced  sprays  calculated  to  have  D^'s  between  20  and  40  microns. 
Because  the  volume  mean  diameters  of  the  liquid  oxygen  sprays  are  only 
slightly  larger  than  the  smallest  droplets  being  produced,  correlation  be¬ 
tween  the  calculated  results  and  standard  distribution  equations  is  almost 
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impossible.  Another  problem  arises  because  the  spruys  calculated  by  the 
coaxial  jet  model  have  both  maximum  and  minimum  droplet  sizes,  features 
that  most  distribution  equations  are  not  capable  of  describing  accurately. 

The  upper  limit  distribution  function  (Ref.  13)  accounts  for  a  maximum 
spray  droplet  diameter  and  was  used  to  correlate  the  spray  data  in  Ref.  3. 
A  dimensionless  function  of  droplet  diameter,  y,  is  assumed  to  be  normally 
distributed;  thus, 
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The  cumulative  volume  fraction  is  then  given  by, 
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Tables  of  the  "probability  integral"  may  be  used  to  evaluate  the  latter 
Integra  1  . 

For  a  spray  distribution  which  fits  the  upper  limit  function,  a  plot  of 

D  (D  -D^  vs  v  on  log-probabi  1  i  t  aper  should  yield  a  straight  line, 
max 

The  values  of  5  and  u  can  be  determined  from  such  a  plot  and,  from  them, 
various  mean  droplet  diameters  can  be  calculated  (Ref.  15). 
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The  calculated  distributions  of  liquid  oxygen  spray  mass  are  shown  plotted 
on  that  basis  in  Fig.  8.  It  is  seen  that  about  70  to  80  percent  of  the 
spray  mass  is  fairly  well  fitted  by  a  straight  line,  and  further,  that  the 
two  sprays  are  very  nearly  identical  on  this  basis.  Here  again,  however, 
the  function  fails  to  describe  the  calculated  accentuation  of  spray  mass 
associated  with  the  smallest  droplet  sizes.  Perhaps  this  is  caused  by  the 
arbitrary  suppression,  in  the  atomization  model,  of  only  the  atomization 
rate  (rather  than  both  rate  and  mean  droplet  sizes)  in  the  first  few  A  x- 
lncrements.  On  the  other  hand,  it  may  be  possible  that  a  lowoi  limit  upper- 
limit  distribution  function  would  be  required  for  adequately  correlating 
sprays  resulting  from  this  kind  of  atomization  when  the  relative  velocity 
is  initially  very  high.  Such  a  distribution  function  is  not  known  at  the 
present  time. 

The  percentages  of  the  liquid  oxygen  flowrate  that  were  calculated  to  have 
been  atomized  in  the  nouburning  region  are  shown  in  Fig.  9;  these  were  ob¬ 
tained  by  subtracting  the  percent  of  liquid  oxygen  remaining  in  the  liquid 
jet  from  100  percent.  The  abrupt  changes  in  the  atomization  rates  at  0.6 
and  0.9  inch  correspond  to  the  cessation  of  atomization  because  any  fur¬ 
ther  atomization  would  produce  droplets  larger  than  the  maximum  permitted 
diameter.  Nearly  imperceptible  increases  in  "percent  atomized"  downstream 
of  those  positions  result  from  continued  evaporation  from  the  surface  of 
the  liquid  oxygen  jets. 


Spray  Vaporization.  AIbo  shown  in  Fig.  9  are  the  calculated  percentages 
of  the  injected  liquid  oxygen  that  were  vaporized.  The  comparatively 
rapid  vaporization  in  the  first  1  to  1-1  2  inches  lrom  the  injector,  par¬ 
ticularly  obvious  under  the  conditions  of  test  417,  was  caused  by  the  ini¬ 
tial  lj  high  velocity  and  concentration  gradients.  In  some  cases,  a  few 
of  the  car  1 lesU-f ormed,  smal lest-diameter  spray  groups  were  calculated  to 
have  been  completely  evaporuted  within  this  portion  of  the  chamber. 
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The  calculated  oxygen  concentrations  at  the  outer  boundaries  of  the  non¬ 
burning  spray  zones  are  shown  in  Fig.  10.  The  effect  of  the  higher  gase¬ 
ous  hydrogen  tenperature  in  test  417  i*  clearly  evidenced  by  the  2  to  3 
percent  higher  oxygen  concentration  throughout  the  nonburning  region. 

The  bumps  in  these  curves  are  loosrly  associated  with  positions  of  com¬ 
plete  vaporization  of  the  outermost  spray  groups,  with  positions  where 
the  spray  zone  boundary  velocity  falls  below  the  droplet  velocities  in 
some  of  the  outermost  spray  groups  and  with  the  completion  of  engulfing 
of  the  hydrogen  stream  into  the  spray  zone. 

Further  insight  into  the  spray  vaporization  phenomena  may  be  gained  by 
studying  Fig.  11.  Here,  at  positions  fairly  close  to  the  injector,  the 
oxygen  concentration  in  the  spray  zone  gases  is  seen  to  be  highest  near 
the  liquid  oxygen  jet  surface.  This  is  because  all  of  the  gases  drawn 
into  that  area  have  passed  through  the  outlying  spray  elements  and  have 
evaporated  part  of  their  liquid.  Additionally,  those  outlying  elements 
are  being  continually  exposed  to  new  hydrogen  entering  the  spray  zone. 
Near  the  injector,  then,  the  evaporated  oxygen  tends  to  migrate  in,  to¬ 
ward  the  liquid  oxygen  jet. 

Further  downstream,  approximately  1  to  l-l/2  inches  from  the  injector, 
the  largest  part  of  the  gas  deceleration  has  occurred,  no  new  spray  is 
formed,  and  little  additional  momentum  is  required  for  accelerating  the 
spray  groups.  The  spray  zone  gases  tend  to  remain  associated  with  a  par¬ 
ticular  spray  group  for  longer  times.  Because  the  heat  of  vaporization 
for  that  liquid  oxygen  that  was  vaporized  had  to  come  from  the  gas  stream, 
the  coldest  gases  are  those  nearest  the  liquid  oxygen  jet.  These  gases 
also  surround  the  largest  droplets.  Vaporization  near  the  liquid  oxygen 
jet  is  thus  very  slow,  while  the  smaller,  outlying  spray  droplets  (sur¬ 
rounded  by  warmer  gases)  vaporize  at  moderate  rates.  The  results  are 
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reflected  in  Fig.  11  as  nearly  constant  oxygen  concentrations  in  the  gases 
near  the  liquid  oxygen  jet  and  continuously  increasing  concentrations  vith 
increasing  radii  and  distances  from  the  injector. 


Velocities  and  Temperatures.  Calculated  velocities  at  the  outer  bound¬ 
aries  of  the  spray  zones  and  the  velocities  of  the  liquid  oxygen  spray 
element  droplets  formed  at  a  number  of  specific  distances  from  the  in¬ 
jector  are  plotted  in  Fig.  12.  The  diameters  of  the  liquid  oxygen  drop¬ 
lets  vhen  they  were  formed  are  indicated  on  the  droplet  velocity  curves. 
Under  the  conditions  of  teat  417  (Fig.  12a),  droplets  of  20.55  microns 
initial  diameter  and  smaller  were  calculated  to  be  completely  vaporized 
within  the  nonburning  region;  the  positions  of  droplet  consumption  are 
indicated  by  simply  terminating  the  velocity  plots.  No  spray  groups  were 
calculated  to  have  been  completely  gasified  in  the  nonburning  region  for 
the  conditions  of  test  431  (Fig.  12b). 

The  calculated  spray  zone  droplet  and  associated  mixed  gas  temperatures 
are  plotted  in  Fig.  13  for  the  same  spray  elements  whose  velocities  were 
plotted  in  Fig.  12a.  To  facilitate  comparisons,  all  curves  are  labeled 
with  the  initial  droplet  diameters  of  the  spray  elements.  It  is  seen 
that  near  the  injector,  as  noted  earlier,  the  gases  cooled  by  oxygen  evap¬ 
orating  into  them  are  carried  radially  inward  toward  the  liquid  oxygen 
jet.  Beyond  1  inch  from  the  injector,  the  gas  temperatures  all  de¬ 
crease  gradually  as  some  of  their  sensible  heat  is  used  to  supply  heat 
of  vaporization.  The  rate  of  gas  temperature  change  is  highest  in  the 
outermost  elements,  i.e.,  those  which  contain  the  smallest,  fastest- 
evaporating  droplets. 

The  liquid  oxygen  droplets  are  calculated  to  be  heated  to  equilibrium 
temperatures  and  undergo  quasisteady-state  evaporation.  Some  of  the 
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largest  droplet  size  groups  had  not  yet  reached  an  equilibrium  tempera¬ 
ture  before  they  entered  the  combustion  region.  The  equilibrium  tem¬ 
peratures  are  seen  to  vary  slightly  with  differences  in  droplet  sizes, 
gas  temperatures,  and  convective  conditions. 

Qualitatively  similar  results  were  calculated  for  conditions  of  test  431. 
The  equilibrium  droplet  temperatures  were  clustered  around  200  to  210  R 
as  opposed  to  the  220  to  230  R  range  shown  in  Fig.  13.  The  160  R  lower 
gas  temperature  was  the  largest  factor  in  causing  this  difference. 


Combustibility.  Examination  of  Fig.  10,  11,  and  12  reveals  that  the  ulti¬ 
mate  satisfaction  of  the  combustibility  conditions  was  controlled  by  dif¬ 
ferent  criteria  under  the  two  different  sets  of  test  conditions.  For 
test  417,  with  high  hydrogen  injection  temperature  and  velocity,  enough 
oxygen  was  evaporated  approximately  2  inches  from  the  injector  to  fulfill 
the  flammability  criterion,  but  the  spray  zone  boundary  gas  velocity  was 
nearly  400  ft/sec.  The  proportionality  constant  in  Eq.  23  was  thus  se¬ 
lected  to  force  to  be  equal  to  the  calculated  spray  zone  boundary  vel¬ 
ocity  at  the  position  of  abrupt  velocity  changes  observed  in  the 
experimental  data  (Ref.  l). 

For  test  431,  just  the  opposite  situation  was  obtained.  The  spray  zone 
boundary  velocity  dropped  below  the  calculated  turbulent  flame  speed  at 
about  2  inches  from  the  injector,  but  the  oxygen  content  of  the  spray 
zone  gases  was  only  about  5  to  6  mole  percent.  The  value  of  the  flam¬ 
mability  limit  was  thus  »et  equal  to  the  spray  zone  boundary  oxygen  con¬ 
centration  at  the  position  of  abrupt  experimental  velocity  change  of  this 
tent. 
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Evaporation  of  a  liquid  droplet  will  cease  when  the  partial  pressure  of 
the  liquid  species  in  the  gas  stream  becomes  equal  to  that  specie^  vapor 
pressure  at  the  droplet  temperature (Eq.  12).  If  oxygen/hydrogen  rockets 
are  operated  with  very  high  chamber  pressures  and/or  with  relatively  low 
propellant  injection  temperatures,  it  may  not  be  possible  to  satisfy  the 
combustibility  conditions.  The  pertinent  relationships  are  shown  in 
Fig.  14,  a  plot  of  the  maximum  achievable  oxygen  mole  fraction  in  gaseous 
hydrogen  as  determined  from  equilibrium  liquid  oxygen  droplet  temperature 
and  chamber  pressure.  The  maximum  equilibrium  droplet  temperature  is 
equal  to  the  gas  temperature,  even  if  the  liquid  were  injected  at  a  higher 
temperature  than  the  gas.  For  the  purposes  of  this  discussion  then,  the 
abscissa  of  Fig.  14  may  be  entered  using  hydrogen  injection  temperature.* 
Coaxial  jet  injection  of  liquid  oxygen  and  gaseous  hydrogen  with  injec¬ 
tion  temperatures  of  183  R  into  a  rocket  operating  at  700  psia  chamber 
pressure,  for  example,  could  not  meet  the  combustibility  requirements 
used  in  the  current  model.  Figure  14  shows  6.5  mole  percent  oxygen  to 
be  the  maximum  attainable.  This  situation  can  be  ameliorated  only  by 
augmenting  the  liquid  oxygen  vaporization  by  supplying  heat  from  some 
other  source  not  considered  in  this  model. 

Two  heat  sources  not  accounted  for  in  the  model  are  radiation  fr^  'he 
chamber  full  of  high-temperature  gases  and  recirculation  of  combustion 
gases.  Radiant  energy  transfer  to  the  nonburning  region  was  examined 
for  the  conditions  of  test  417  in  a  simplified  manner;  it  was  calculated 
that  as  much  as  1  percent  of  the  liquid  oxygen  injection  rate  could  be 

*Figure  13  suggests  that  a  value  of  gas  tesq>erature  10  to  20  R  lower  than 
hydrogen  injection  temperature  would  be  more  realistic.  This  tempera¬ 
ture  drop  provides  the  sensible  heat  of  the  warming  droplets,  latent  heat 
of  vaporization,  and  superheat  for  the  vapors. 
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evaporated  by  radiation.  At  an  injection  mixture  ratio  of  5»  however, 
the  oxygen  content  of  the  spray  zone  gases  would  be  raised  by  an  average 
of  only  0.3  mole  percent.  Under  other  operating  conditions,  this  effect 
could  conceivably  have  more  influence  and  needs  to  be  taken  into  account. 

The  likelihood  of  combustion  gas  recirculating  was  indicated  by  the  value 
of  the  Craya-Curtet  number,  previously  described.  The  calculated  values 
of  that  variable  for  the  two  tests  analyzed  are  shown  in  Fig.  15*  In 
both  cases,  the  Craya-Curtet  number  exceeded  the  critical  value  of  3A 
well  before  the  combustibility  conditions  had  been  satisfied;  the  omis¬ 
sion  of  combustion  gas  recirculation  from  the  model  thus  appears  to  be 
a  valid  simplification  for  these  cases. 

Chamber  Pressure.  The  majority  of  decrease  in  gas  stream  momentum  as  the 
gases  are  decelerated  is  absorbed  in  accelerating  the  liquid  oxygen  spray. 
The  remainder  results  in  increasing  chamber  pressure  along  the  axial 
length.  The  model  calculations  showed  a  rapid  rise  in  chamber  pressure 
in  the  first  l/ 2  inch  of  the  chamber  and  a  gradual  rise  in  pressure 
throughout  the  remainder  of  the  nonburning  region.  Certainly  some  pres¬ 
sure  rise  should  be  expected,  but  in  the  absence  of  direct  experimental 
evidence,  it  is  not  possible  to  state  how  realistic  the  model  predictions 
are  in  this  regard.  Possibly,  the  model  overestimates  the  amount  of  gas 
needed  to  be  decelerated  for  providing  spray  momentum,  and  thus  overempha¬ 
sizes  the  pressure  rise  in  this  region.  This  subject  will  be  discussed 
further  in  the  next  section. 
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Combustion  Region 

The  r-^nburning  regions  vere  replaced  by  combustion  regions  at  3.25  and 
3.20  inches  from  the  injector  in  aodel  calculations  based  on  the  condi¬ 
tions  of  tests  417  and  431 ,  respectively. 


Structure .  The  combustion  region  model  resulted  in  structural  character¬ 
istics  shown  in  Fig.  16  for  test  417.  This  figure  is  a  direct  follov-on 
to  Fig.  4.  Comparison  with  Fig.  1  shows  the  calculated  structure  to  be 
roughly  similar  to  that  initially  hypothesised;  the  calculated  spray  zone 
spreading  is  more  abrupt  and  the  nonspray-bearing  gas  mantle  is  much  less 
prominent  than  anticipated. 

The  abrupt,  rapid  spreading  of  the  propellant  streams  to  completely  fill 
the  elemental  chamber  flow  channel  is  an  idealization  required  by  the 
stepwise,  downstream-progression  model  employed.  In  an  actual  subsonic 
flow,  the  propellant  flow  structure  upstream  of  the  flame  front  should  be 
affected  by  the  impending  combustion  region  expansion  and  the  spreading 
made  less  abrupt  than  that  calculated. 

Intuitively,  the  structure  shown  in  Fig.  16  suggests  that  the  choice  of 
this  combustion  region  model  based  on  a  well-mixed,  one-dimensional  spray 
combustion  zone  being  surrounded  by  a  nonspray-bearing  gas  mantle  is 
probably  much  less  realistic  than  one  based  on  a  very  oxygen-rich,  dense 
spray  core  being  surrounded  by  combustion  gases  having  decreasing  spray 
population  with  increasing  distance  from  the  axial  centerline.  Having 
selected  the  former  basis  for  the  model  formulation,  however,  other  cal¬ 
culated  results  will  be  examined  before  critical  conclusions  are  drawn. 
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Combust  ion  Gas  Velocity.  Inuediately  upon  passing  through  the  flame 
front  at  the  leading  edge  of  the  combustion  region,  the  spray  zone  gases 
were  reacted.  An  abrupt  increase  of  volumetric  flowrate  resulted;  this, 
in  turn,  caused  radial  expansion  until  the  propellant/combustion  gas  flow 
completely  filled  the  flow  channel.  Then,  higher  axial  gas  velocities  were 
attended  by  decreased  chamber  pressure.  The  existing  liquid  oxygen  drop- 
lets'  vaporization  rates  were  greatly  enhanced  as  a  result  of  both  higher 
gas  temperatures  and  greater  convection.  This  further  raised  the  combus¬ 
tion  gas  velocity.  The  gas  velocity  soon  became  high  enough  so  that  liq¬ 
uid  oxygen  jet  atomization  by  the  surface  shear  mechanism  vas  resumed. 
These  processes  are  thus  seen  to  be  interrelated  in  a  manner  that  causes 
escalation  or  bootstrapping  to  quite  high  combustion  gas  velocities. 


Comparison  Vith  Experimental  Velocity  Data.  The  experimentally  derived 
velocity  and  pressure  data  presented  in  Ref.  1  for  tests  417  and  431  are 
reproduced  in  Fig.  17a  and  17b,  respectively.  The  calculated  combustion 
gas  velocities  (and  other  parameters  which  are  discussed  below)  are  also 
plotted  to  facilitate  direct  comparisons.  It  is  seen  that  the  calculated 
combustion  gas  velocity  jumps  agree  quite  veil  vith  those  observed  experi¬ 
mentally.  This  suggests  that  the  early  portions  of  the  combustion  model 
adequately  describe  the  initial  bootstrapping  processes. 


Liquid  Oxygen  Jet  Atomization.  Under  the  conditions  of  test  417,  the 
atomization  process  vas  calculated  to  be  resumed  immediately  downstream 
of  tiie  flame  front.  The  intially  high  combustion  gas  velocity  resulted 
immediately  m  quite  small  liquid  oxygen  droplets;  less  than  l/4  inch  in¬ 
to  the  combustion  region,  the  produced  spray  was  calculated  to  be  wholly 
vaporized  and  reacted  within  the  same  Ax  increment  that  it  vas  atomized. 
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The  residual  liquid  oxygen  jet,  containing  45  percent  of  the  liquid  oxy¬ 
gen  injection  rate,  was  thus  found  to  be  completely  atomized  in  only  l/2 
inch. 

Because  of  the  lover  hydrogen  injection  rate  and  the  lover  proportion  of 
oxygen  evaporated  in  the  nonburning  region  of  test  431,  a  lesser  jump  in 
gas  velocity  was  calculated  to  occur  at  the  beginning  of  the  combustion 
region  than  was  found  for  test  417.  The  resisaption  of  liquid  oxygen  jet 
atomization  was  delayed,  therefore,  until  accelerated  droplet  vaporiza¬ 
tion  increased  the  combustion  gas  velocity  even  more.  Droplet  sizes  near 
the  maximum  were  thus  produced  at  first. 

Continued  vaporization  led  to  a  gradual  reduction  in  the  produced  droplet 
sizes,  but  it  was  nearly  1  inch  into  the  combustion  region  before  the 
atomization  rate  became  comparable  with  that  found  in  the  combustion  re¬ 
gion  for  test  417.  These  relationships  are  seen  in  Fig.  18.  The  calcu¬ 
lated  ends  of  the  residual  liquid  oxygen  jets  correspond  veil  with  the 
experimentally  observed  distances  to  which  velocity  as  low  as  liquid  oxy¬ 
gen  injection  velocity  penetrated  in  Fig.  17. 

Two  representative  examples  of  the  droplet  weight  distributions  are  shown 
in  Fig.  19.  The  first,  Fig.  19a, is  for  the  sprays  traversing  the  flame 
fronts  the  second,  Fig.  19b,  is  for  the  sprays  remaining  just  after  com¬ 
pletions  of  liquid  oxygen  jet  atomization.  Comparing  Fig.  19a  with  Fig. 
7b  reflects  the  complete  evaporation  of  several  liquid  oxygen  spray  group 
sizes  in  the  nonburning  region  for  test  417,  so  that  D^q  increased  and 
the  distribution  favored  the  larger  drop  sizes,  but  there  was  only  par¬ 
tial  vaporization  of  all  spray  group  sizes  for  test  431,  so  that  D^q 
simply  decreased.  These  changes  do  not  appear  to  be  predictable,  a 
priori,  from  the  injection  conditions. 
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Comparing  Fig.  19a  and  19b,  it  is  seen  that  the  spray  throughout  the  com¬ 
bustion  region  for  test  417  can  be  described  by  the  Nukiyama-Tanasawa  dis¬ 
tribution  function.  Because  the  new  spray  atomized  from  the  jet  was  found 
to  evaporate  as  soon  as  it  was  formed,  all  the  spray  shown  for  test  417  in 
these  figures  had  its  origin  in  the  nonburning  region.  Ingebo  (Ref.  16) 
has  observed  experimentally  that  a  burning  ethyl  olcohol  spray  similarly 
fits  the  Nukiyama-Tanasava  function  at  several  distances  from  a  rocket 
injector.  The  liquid  oxygen  spray  in  the  combustion  region  for  test  431 
was  not  nearly  so  simply  behaved:  nearly  83  percent  of  its  spray  shown 
in  Fig.  19b  originated  in  •.  e  combustion  region.  Again,  neither  the  dis¬ 
tributions  nor  mean  droplet  sizes  of  these  sprays  appear  to  be  predictable 
by  any  means  other  than  this  model. 


Spray  Droplet  Behavior.  Calculated  velocities  for  some  droplet  size  groups 
are  plotted  in  Fig.  17a  and  17b.  The  velocities  were  comparable  with  those 
calculated  in  Appendix  B  and  were  well  behaved.  The  largest  formed  drop¬ 
let's  velocities  described  a  lower  boundary  for  the  experimental  velocity 
data  points. 

Droplet  equilibrium  temperatures  of  240  R  were  calculated  for  the  condi¬ 
tions  of  both  tests  in  the  early  portions  of  the  combustion  regions.  These 
gradually  dropped  to  235  R  at  the  beginning  of  nozzle  convergence,  in  di¬ 
rect  agreement  with  the  behavior  predicted  in  Appendix  B,  Fig.  B-l  and  B-2. 


Nonspray-Bearing  Gas  Mantle  Behavior.  The  gaseous  flow  diverted  out  of 
the  nonburning  liquid  oxygen  spray  zone  and  around  the  flame  front  con¬ 
tained  enough  oxygen  to  be  flanmable  under  the  conditions  of  test  417, 
but  not  under  the  conditions  of  test  431.  The  model  calculations  thus 
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ascribed  reacted,  combustion-gas  properties  to  this  gas  stream  in  the 
former  case  and  unreacted,  mixed-gas  properties  in  the  latter  case.  The 
calculated  gas  stream  velocities  were  correspondingly  high  for  test  417 
and  low  for  test  431,  as  seen  in  Fig.  17.  Because  there  vere  no  experi¬ 
mental  velocity  data  points  as  high  as  the  calculated  gas  stream  velocity 
for  test  417,  it  must  be  concluded  that  either  the  radiation  from  that 
thin,  reacted  mantle  was  too  weak  to  affect  the  streak  film  exposures  or 
the  mantle  is  fictitious. 

Mixing  of  the  model-calculated  gas  stream  flowrate  into  the  spray  combus¬ 
tion  zone  at  a  linear  rate  with  increasing  distance  from  the  injector, 
however,  resulted  in  a  combustion  gas  velocity  profile  that  modeled  the 
highest  experimental  velocity  data  very  veil.  This  is  seen  in  the  cham¬ 
ber  length  between  3  inches  from  the  injector  and  the  beginning  of  the 
nozzle  convergence  (Fig.  17). 

/ 

To  determine  whether  the  gas  velocity  profile  might  be  determined  by 
unevaporated  oxygen  spray  distributed  along  the  chamber  (Fig.  18),  model 
calculations  were  made  with  no  spray  zone  gases  diverted  into  a  gas 
stream  mantle.  The  results  (Fig.  20)  suggest  that  the  observed  gradual 
approach  to  the  nozzle  inlet  combustion  gas  velocity  is  almost  wholly 
attributable  to  turbulent  mixing  of  unlike  gas  streams.  As  noted  before, 
the  actual  situation  might  require  mixing  a  strongly  oxidizer-rich  core 
radially  outward  rather  than  a  fuel-rich  mantle  radially  inward.  The 
overall  effect  would  be  the  same;  (her  a  sizeable  fraction  of  the  chamber 
length,  the  propellant  combustion  is  controlled  by  turbulent  mixing 
processes. 

This  result  of  the  present  model  calculations  is  in  direct  agreement  with 
the  results  of  related  experimental  and  analytical  studies  at  NASA,  Lewis 
Research  Center  (Ref.  18  and  19). 
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Chamber  Pressure.  The  calculated  chamber  pressures  are  seen  to  decrease 
rapidly  immediately  downstream  of  the  flame  fronts.  These  drops  result 
from  continuity  of  momentum.  The  increasing  momenta  of  accelerating 
gases  and  spray  can  come  only  from  a  decreasing  chamber  pressure. 

Following  that  precipitous  drop,  the  calculated  chamber  pressures  agree 
well  with  the  measured  pressures.  Too  few  pressures  were  measured  in  the 
nonburning  region  and  early  combustion  region  to  define  the  actual  pres¬ 
sure  profile,  but  those  obtained  suggest  smoother  pressure  profiles  than 
those  calculated.  Perhaps  the  apparent  disagreement  stems  from  the  ex¬ 
perimental  measurements,  e.g.,  the  pressure  tap  geometry  or  flow  patterns 
may  have  affected  the  data  and  (because  the  data  are  time  and  location 
averaged)  fluctuations  in  individual  flame  front  position  and  variations 
in  those  positions  among  all  the  jets  would  tend  to  smooth  and  flatten 
the  pressure  profile.  On  the  other  hand,  the  model  could  be  forced  to 
give  a  more  nearly  constant  pressure  profile  by  assuming  more  drastic 
suppression  of  atomization  near  the  injector,  considerably  higher  turbu¬ 
lent  flame  speeds  and  different  flammability  limits  than  those  used.  The 
calculated  gas  velocities  would  no  longer  match  the  experimental  streak 
velocities  upstream  of  the  flame  front,  however.  More  quantitative  know¬ 
ledge  in  this  regard  will  probably  be  dependent  upon  both  refined  pres¬ 
sure  measurements  and  determination  of  f  1  amiability  limits  and  flame 
velocities  under  appropriate  rocket  conditions. 

Concluiuons  From  the  Best-Fit  Cases 

The  integrated  atomization  and  spray  behavior  formulation  permits  des¬ 
cription  of  the  entire  propellant  flow  structure  from  specifications  of 
the  injector  and  chamber  configurations  and  propellaut  injection  conditions. 
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The  potential  utility  of  the  combustion  model  is  thus  greatly  enhanced. 
Further,  the  formulation  should  be  directly  applicable  to  studies  of  simi¬ 
lar  atomization  processes  with  other  fluids,  whether  or  not  evaporation 
is  significant.  For  extension  to  other  rocket  combustors,  the  formula¬ 
tion  may  need  to  be  expanded  to  include  nonreacting  gas  circulation,  com¬ 
bustion  gas  recirculation,  and  continued  atomization  by  other  processes, 
such  as  jet  instability  and  centrifugal  forces  imposed  by  svirlers. 

The  model  calculations  offer  definitive,  convincing  evidence  for  the  exis¬ 
tence  of  the  hypothesized  plane  flame  front  standing  at  a  fixed  (or  quasi- 
fixed)  position  in  the  liquid  oxygen  spray  zone.  Under  the  imposed  model 
conditions,  any  penetration  of  combustion  gases  or  flame  into  this  zone 
at  some  upstream  position  would  either  be  swept  downstream  (by  the  gas 
flow  velocity  exceeding  the  local  turbulent  flame  speed)  or  would  rapidly 
escalate  through  increased  spray  vaporization,  chamber  filling,  gas  ac¬ 
celeration  and  enhanced  liquid  oxygen  jet  atomization  to  the  some  pheno¬ 
mena  described  by  the  model  but  at  the  upstream  position.  This  seemingly 
adamant  interpretation  is  dictated  by  the  model  conditions  of  relatively 
low  contraction  ratio  (which  permits  moderately  high  combustion  gas  veloc¬ 
ities  even  though  only  fractional  amounts  of  the  liquid  oxygen  have  been 
vaporized,  and  leads  to  combustion  reaction  and  atomization  mutually  aug¬ 
menting  each  other)  and  unappreciable  combustion  gas  recirculation. 

This  interpretation  is  not  believed  to  be  incompatible  with  nor  to  refute 
other  investigators'  interpretations  of  their  experimental  observations 
in  which  it  is  suggested  that  the  atomizing  liquid  oxygen  jet  and  liquid 
oxygen  spray  zone  are  intimately  surrounded  by  reacting  combustion  gases 
all  the  way  upstream  to  the  injector  (Ref.  19  and  21 ).  Their  experiments 
were  carried  out  in  relatively  high-contraction-ratio  combustors  and 
under  conditiors  conducive  to  combustion  gas  recirculation,  so  that  the 
bootstrapping  phenomena  should  not  have  occurred. 
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The  combustion  region  model  appears  to  constitute  an  adequate  although 
somevhat  simplified  representation  of  the  spray  gasification  and  trans¬ 
verse  mixing  processes.  The  one-dimensional  model  of  the  nozzle  flow 
processes  appears  to  describe  them  accurately  to  the  nozzle  throat  when 
a  nozzle  erosion  factor  is  applied  to  account  for  the  cross-sectional 
area  variation  resulting  from  pyrex  melting. 

MODEL  PREDICTIONS  FOR  OTHER  MA INSTAGE  CONDITIONS 


Model  calculations  were  carried  out  in  which  controlled  variations  were 
made  in  the  initial  conditions,  while  holding  the  values  of  all  arbitrarily 
evaluated  variables  equal  to  those  used  in  the  best-fit  cases.  Chamber 
pressure,  mixture  ratio,  hydrogen  injection  temperature,  and  liquid  oxygen 
injection  temperature  were  individually  varied  from  those  listed  in  Table  1 
for  test  431.  The  major  comparison  bases  for  evaluating  the  predicted  ef¬ 
fects  were  the  location  in  the  chamber  where  the  combustibility  criteria 
were  satisfied,  and  which  of  those  criteria  was  the  last  to  be  fulfilled. 

Raising  the  chamber  pressure  by  increasing  both  propellant  flowrates  at 
constant  injection  mixture  ratio  resulted  in  the  calculated  position  of 
the  flame  front  being  closer  to  the  injector  face,  i.e.,  3.2  inches  at 
435  psia,  2.8  inches  at  300  psia,  2.3  inches  at  700  psia,  and  1.7  inches 
at  800  psia.  At  the  lower  pressures,  this  effect  resulted  from  somewhat 
finer  atomization  because  of  the  increased  gas  density  at  higher  pressures. 
At  800  psia,  supercritical  heating  of  the  smallest  liquid  oxygen  droplets 
was  encountered.  Their  abrupt  gasification  caused  the  combuatibility 
criteria  to  be  satisfied  appreciably  earlier  than  they  were  at  700  psia. 
Hydrogen  injection  temperature  was  high  enough  to  avoid  altogether  the 
vaporization  difficulties  noted  earlier  in  relation  to  Fig.  14.  Calculated 
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boundary  velocities  dropped  below  the  turbulent  flame  speeds  earlier  than 
the  flamnable  oxygen  concentrations  were  reached  in  all  of  these  cases. 

The  Craya-Curtet  number  at  the  flame  front  was  less  than  its  critical 
value  in  the  800  psia  case,  implying  that  combustion  gas  recirculation 
could  invalidate  the  model  predictions  at  that  and  higher  pressures. 

The  injection  mixture  ratio  was  varied  by  raising  the  hydrogen  flowrate 
and  lowering  the  liquid  oxygen  flowrates  by  amounts  that  would  give  total 
flowrate  satisfying  100  percent  c*  efficiency  at  a  fixed  chamber  pressure. 
The  flame  front  position  moved  downstream  with  increasing  mixture  ratio 
and  was  c  ntrolled  by  the  attainment  of  flammable  mixtures.  The  flame 
front  moved  upstream  with  decreasing  mixture  ratio,  in  general,  but  at 
mixture  ratios  of  approximately  5.  flanmability  could  be  achieved  only 
at  sizable  distances  (10  to  12  inches)  downstream  of  the  injector.  The 
reason  for  this  anomalous  behavior  was  not  discerned.  Three  position 
values  were:  (l)  2.5  inches  at  MR  =  3.  (2)  3.2  inches  at  MR  =  7,  and 
(3)  3.9  inches  at  MR  =  12.  Combustibility  was  controlled  at  that  lowest 
mixtu.e  ratio  by  the  velocity  criterion. 

Liquid  oxygen  injection  temperature  was  varied  from  160  to  210  R  and  the 
calculated  positiou  of  the  flame  front  moved  steadily  upstream  as  the  tem¬ 
perature  was  increased.  Some  values  were:  approximately  6  inches  at 
160  R,  3.4  inches  at  173  to  180  R,  2.4  inches  at  190  R,  and  1.9  inches 
at  210  R.  Reaching  flammable  oxygen  concentrations  controlled  combusti¬ 
bility  at  temperatures  below  200  R,  while  deceleration  to  turbulent  flame 
velocity  controlled  those  cases  having  200  R  or  warmer  liquid  oxygen. 

The  flame  front  position  similarly  moved  upstream  with  hydrogen  injection 
temperature  increases.  At  400  R,  the  combustibility  was  controlled  by 
the  velocity  criteria,  which  was  satisfied  at  2.7  inches.  The  computer 
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model  encountered  computational  difficulties  in  the  nonburning  spray  zone 
with  TOO  R  hydrogen  und  would  not  converge  on  a  solution.  With  hydrogen 
colder  than  that  reported  for  test  431 »  the  flame  front  moved  downstream. 
While  200  R  hydrogen  should  eventually  produce  a  combustible  mixture 
(Fig.  14),  the  liquid  oxygen  evaporation  was  so  slow  that  it  is  question¬ 
able  whether  combustibility  could  be  achieved  in  the  available  chamber 
length;  the  highest  spray  zone  oxygen  concentration  at  5  inches  from  the 
injector  was  only  about  3  percent. 


Comparison  With  Experimental  Data 

Several  other  tests  were  reported  in  Ref.  1  to  have  been  conducted  with 
the  coaxial  jet  injector.  While  the  combustion  field  velocity  data  were 
incomplete  and  a  flame  front  position  cannot  be  ascribed  to  each  of  those 
tests,  comparisons  can  be  based  on  the  lengths  of  the  liquid  oxygen  jet/ 
spray  visible  in  the  high-speed  Fastax  motion  pictures.  The  experimental 
condition  for  all  the  tests,  together  with  estimates  of  the  liquid  oxygen 
jet  lengths,  are  sunmarized  in  Table  2. 

Chamber  pressure  variation  among  the  tests  was  quite  modest  and  most  likely 
did  not  influence  the  observed  liquid  oxygen  jet  lengths  appreciably. 

Liquid  oxygen  injection  temperatures  were  not  measured  and  so  cannot  be 
discussed.  The  combustor  starting  procedure  was  such,  however,  that  the 
steady -state  mains  cage  liquid  oxygen  temperatures  should  have  been  com¬ 
parable  among  the  tests. 

Regarding  the  other  parametric  variations,  the  experiments  do  not  appear 
to  support  the  predicted  anomalous  behavior  at  mixture  ratios  in  the 
vicinity  of  5»  nor  the  predicted  difficulty  in  achieving  combustible 
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mixtures  vith  200  R  hydrogen.  There  may,  of  course,  be  an  interplay  bt- 
tveen  the  tvo  effects  masked  by  the  particular  combinations  of  experimental 
temperatures  and  mixture  ratios  encountered.  However,  in  view  of  the  pre¬ 
dicted  sensitivity  of  the  flame  front  position  to  liquid  oxygen  injection 
temperature,  a  few  degrees  difference  in  that  variable  among  the  tests  may 
have  exhibited  an  unexpectedly  large  influence  on  the  results.  For  that 
reason,  it  was  thought  impossible  to  examine  the  existing  experimental 
data  satisfactorily  with  a  comprehensive  parametric  model  study. 


MODEL  HIEDICTIONS  FOR  H1ESTAGE  CONDITIONS 

It  was  reported  in  Ref.  1  that  the  liquid  oxygen  jet/sprays  were  visible 
during  prestage  for  only  about  l/2  to  3/4  inch  from  the  injector,  and  that 
strong  combustion  gas  recirculation  was  evident  in  the  motion  pictures. 
Because  the  prestage  start  transient  lasted  for  about  0.4  to  0.3  second 
and  approached  steady-state  conditions  near  the  end  of  that  time,  model 
calculations  were  programed  for  two  sets  of  prestage  conditions  to  see 
if  their  predictions  would  correspond  with  the  experimental  observations. 

Altering  the  propellant  flowrates  from  those  of  test  431  to  obtain  an  in¬ 
jection  mixture  ratio  of  3  at  260  psia  chamber  pressure  resulted  in  a  pre¬ 
dicted  flame  front  at  3.9  inches  from  the  injector  without  recirculation 
of  combustion  gases.  Similarly,  an  oxygen  flowrate  increase  to  obtain  a 
mixture  ratio  of  3  at  363  psia  chamber  pressure  (corresponding  to  the  mid¬ 
point  between  prestage  and  mainstage)  gave  a  predicted  flame  front  at 
3.05  inches  without  recirculation.  In  both  cases,  deceleration  to  turbu¬ 
lent  flame  velocity  controlled  the  positions  of  the  combustion  fronts. 
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These  results  are  diametrically  opposed  to  the  experimental  observations. 
They  apparently  can  only  be  explained  by  assuming  that  one  or  both  of  the 
propellant  streams  were  actually  much  warmer  than  assumed  in  the  model 
calculations  so  that  the  liquid  oxygen  was  atomized,  accelerated,  and  va¬ 
porized  considerably  faster  than  calculated.  Again,  it  is  unfortunate 
that  liquid  oxygen  manifold  temperatures  were  not  measured.  Also,  fast- 
response  instrumentation  would  have  been  required. 

PREDICTED  OSCILLATORY  COMBUSTION 

A  possibility  that  the  flame  front  might  fluctuate^  rather  than  stand  at  a 
fixed  position  was  suggested  by  model  calculations  for  a  combustor  having 
a  different  shape  and  contraction  ratio  than  the  transparent  two-dimensional 
combustor.  Results  with  200  R  hydrogen  predicted  a  flame  front  standing 
several  inches  from  the  injector,  but  the  expans ion-deccleration-chamber 
area  relationships  were  such  that  strong  combustion  gas  recirculation  was 
indicated.  The  model  prediction  with  300  R  hydrogen  was  that  the  flame 
front  would  be  appreciably  closer  to  the  injector,  but  recirculation  would 
not  be  likely.  Thus,  recirculation  of  enough  combustion  gas  to  warm  the 
hydrogen  by  100  R  would  not  only  permit  the  flame  front  to  move  upstream 
but  would  simultaneously  disrupt  the  combustion  gas  recirculation.  The  con¬ 
tinued  injection  of  fresh,  cold  hydrogen  would  then  push  the  flame  front 
back  ('ownstream  and  reopen  the  recirculation  paths.  Such  an  oscillating 
combustion  process  has  not  been  observed  in  current  combustors  of  conven¬ 
tional  design. 
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CONCLUSION 

It  has  been  found  possible,  using  conceptual  hypotheses  derived  from 
transparent  model  rocket  combustor  experiments,  to  construct  a  semiempiri- 
cal,  analytical  model  of  the  steady-state  combustion  processes  following 
coaxial  jet  injection  of  liquid  oxygen  and  gaseous  hydrogen.  Further,  by 
judicious  selection  of  assumptions  and  particular  arbitrary  constants,  it 
has  been  possible  to  obtain  substantial  agreement  between  the  model's 
predictions  and  experimental  results  from  two  of  the  experiments.  The 
other  experimental  results  have  not  fully  corroborated  the  accuracy  of 
analytically  predicted  trends,  but  insufficient  data  were  obtained  to  es¬ 
tablish  or  deny  the  model's  validity. 

Several  aspects  of  the  developed  rocket  combustion  model  are  unique.  Per¬ 
haps  the  most  striking  feature  is  the  interpretation  of  the  experimental 
results  in  terms  of  a  standing,  plane  flame  front.  This  concept  has  been 
determined  conclusively  to  be  inappropriate  for  systems  in  which  both  pro¬ 
pellants  are  injected  as  liquids,  but  appears  from  this  work  to  be  valid 
for  liquid-gas  combinations  if  combustion  gas  recirculation  is  insignificant. 

The  use  of  fundamental  combustion  parameters  (f lanability  limits  and  tur¬ 
bulent  flame  speeds)  in  rocket  combustion  models  is  unusual;  they  are  here 
related  directly  to  the  flame-front  phenomenon.  If  continued  uses  are  to 
be  made  of  coaxial  jet  injectors  under  conditions  similar  to  those  studied 
here,  it  would  prove  very  helpful  to  analytical  model  studies  to  have  more 
accurate  knowledge  of  these  parameters  under  appropriate  rocket  pressure, 
temperature,  and  flow  conditions. 
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Another  unique  feature  of  the  current  model  is  the  quantitative  inclusion 
of  the  atomization  process.  This  is  the  first  time  to  our  knowledge  that 
this  has  been  attempted.  The  significance  of  this  accomplishment  was  some¬ 
what  reduced  by  the  necessity  of  using  arbitrary  constants.  Currently,  the 
nonburning  region  model  is  being  programed  to  calculate  liquid  wax  atomi¬ 
zation  by  an  airstream  in  anticipation  that  a  single  set  of  constants  can 
be  found  to  reproduce  the  variety  of  atomization  data  in  Ref.  3>  md  that 
those  can  be  related  to  the  constants  used  here  for  liquid  oxygen  gaseous 
hydrogen . 

One  reason  for  developing  a  steady-state  combustion  model  is  to  determine 
if  combustion  efficiency  will  be  degraded  by  large  propellant  droplets 
passing  through  the  nozzle  without  being  evaporated.  By  including  the 
atomization  process  in  the  model,  direct  knowledge  is  gained  of  what  the 
largest  droplet  sizes  are,  how  much  of  the  spray  masB  is  associated  with 
them,  and  equally  important,  where  they  are  formed.  Thus,  if  undesirably 

coarse  spray  is  encountered,  a  much  more  enlightened  injector  modifica¬ 
tion  program  can  be  undertaken  than  if  only  the  undesirable  droplet  sizes 
were  known. 

Turbulent  mixing  processes  between  unlike  gas  streams  in  the  combustion 
region,  although  treated  in  an  approximate  manner,  have  been  shown  to 
limit  the  propellant  combustion  rates  through  much  of  the  chamber  length. 
The  selected  turbulent  mixing  of  a  gas  mantle  into  a  well-mixed  spray  com¬ 
bustion  zone  is  possibly  less  factual  than  is  the  mixing  of  an  oxygen- 
rich,  dense-spray  core  into  a  surrounding  fuel-rich,  dilute-spray  flow. 
Further,  a  more  rigorous  quantitative  treatment  of  turbulent  mixing  would 
be  required  for  studying  effects  on  combustion  efficiency  of  element  spac¬ 
ing  or  injection  density  changes. 
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Radiant  energy  transfer  and  combustion  gas  recirculation,  neglected  in 
this  model,  appear  to  be  important,  if  not  controlling  processes,  with 
some  injector  and  chamber  geometries  and  under  some  operating  conditions. 
While  these  phenomena  might  profitably  be  included  in  future  model  devel¬ 
opments,  the  fact  that  such  a  statement  can  be  made  demonstrates  one  area 
of  usefulness  of  the  present  model.  Comparison  between  analytical  pre¬ 
dictions  and  experimental  observations  with  a  variety  of  systems  can  in¬ 
crease  greatly  the  current  understanding  of  the  combustion  piocesses. 

Such  qualitative  contributions  are  of  value  through  indirect  influence 
on  problem  solution  across  a  broad  range  of  problems. 

Another  area  of  combustion  model  utility  is  in  predicting  the  effects  of 
making  parametric  changes  in  a  particular  combustor  system.  Quantitative 
effects  can  be  observed.  Of  course,  extension  to  conditions  very  much 
different  from  those  previously  found  to  agree  with  the  model  will  be  sub¬ 
ject  to  some  doubt  unless  limited  experimental  confirmation  is  undertaken. 
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APPENDIX  A 


PHYSICAL,  THERMODYNAMIC,  AND  TRANSPORT  PROPERTIES  FOR  THE 

NONBURNING  REGION 


Data  for  selected  properties  of  hydrogen  and  oxygen  were  curve  fitted  to 
obtain  relationships  for  use  in  the  combustion  model  solutions.  Unless 
a  specific  indication  is  made  of  pressure  or  pressure  dependence,  the 
equations  were  derived  for  properties  at  500  psia  and  the  properties 
were  assumed  not  to  vary  appreciably  with  pressure  within  the  ranges  of 
pressure  encountered  in  the  model  calculations. 


PROPERTIES  OF  LIQUID  OXYGEN 


Density 


pL0X  =  39  88  +  0A19  V  °-001733  T \ 
(Tt<275R) 


(A-l) 


Surface  Tension 

aL0X  =  2,485  X  10~3  "  0,975  X  10"5  Tl 
( 126R  <T^<  160R) 


(A-2) 
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Viscosity 

—  ■  ■  ,, 


The  viscosity  at  zero  pressure  was  arbitrarily  multiplied  by  two  as 
an  approximation  of  the  effect  of  pressure. 


^LOX  "  2  6XP  L  "  11,13  ] 

( 140R  <  <  220R) 


(A-3) 


Specific  Heat  (Constant  Pressure) 


(c  )lnv  =  0.3726  +  2.0482  x  10"^  T, 
p  LOX  /, 


(A-4) 


Vapor  Pressure 


[n.%i  -  ^ 

( IbOR  <  tl  <  278. hR) 


(A-5) 


Heat  of  Vaporization 

The  \alues  at  saturation  pressures  were  curve-fitted; 


H 


v,0o 


128.23 


278.6, 


0.38 


(A-6) 
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GASEOUS  OXYGEN  PROPERTIES 


Specific  Heat  (Constant  Pressure) 


^  Vg<>2 


1.16  -  0.00275  T  ,  (T  <  300  R) 

0.2134  +  (r-:"23r:9)>  (3°°r  <  t  <  740R)  (a-7) 


^Cp^gOf 


0.2489  +  0.964  x 


IQ”5  T  - 


1.555  x  10  ,  (740R  <  T) 


Diffuaivity  into  Hydrogen 

An  expression  recommended  in  Ref.  22  was  used  for  extrapolating  from 
n  value  of  WP  at  273  K  and  1  atmosphere: 

1  A. 


1.4083  x  10~b  T5  2 
P  (T  +  192.6) 


(A-8) 


GASEOUS  HYDROGEN  H10PERTIES 


Viscosity 


*  1.213  x  10-3  +  9.55  x  10"9  T 

(70R  <  T  <  500R) 


(A-9) 
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The  expressions  for  the  following  three  parameters  should  not  be  us~d 
at  temperatures  lover  than  those  indicated. 


Thermal  Conductivity 


kjj  =  0.019  +  1.675  x  10'4  T 
(170R  <  T  <  500R) 


(A-10) 


Specific  Heat  (Constant  Pressure) 


(c  )  „  =  2.095  +  5.08  x  10“3  T  -  4.95  x  10"u  T*  (A-ll) 

r  K 


.-6  „2 


Pg 


(150R  <  T  <  500R) 


Ratio  of  Specific  Heat9 


T  +  2.7 

0.8012T  -  44.4 


/  1  +  2.7  ,  .  v  ( 100  -  P.  1*^2) 

v  0.8012T  -  44.4  "  *  '  v  198  ' 


(T  >  160R) 


( A— 12 ) 
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GASEOUS  MIXTURES  OF  HYDROGEN  AND  OXYGEN 


Mole  Fraction  Oxygen 


0, 


+  15.873 


15873  w  -  14  873  w 

8 


(A- 13) 


Molecular  Weight 


15.873  wg  - 


14.873  w 

2 


(A-H) 


Heat  Capacity  (at  Const  mt  Pressure) 


52  \  (ep)oa  *  2-016  <l  -  (cpV,  ,  4 

(c  )  «  - 5 -  ( A- 1 5 ) 

v  p'g  M 


g 


Viscosity 


The  viscosities  of  mixtures  containing  3i  5»  and  7  mole  percent  oxygen 
were  calculated  by  a  method  recommended  in  Ref.  22  and  fitted  with  the 
following : 


(■ 


1  +  8.5 


(A-  lb) 
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Thermal  Conductivity 

A  correlation  method  recommended  in  Ref.  22  for  predicting  the  thermal 
conductivity  of  gaseous  mixtures  was  used  to  estimate  values  for  3,  5, 
and  7  mole  percent  oxygen  These  were  curve-fitted  by, 

k  =  lv  (1-6.98X  +  32.5  xj!  )  (A-17) 


PROPERTIES  OF  GAS  FILM  AROUND  A  EROPLET 


The  appropriate  preceding  properties  were  evaluated  in  the  gas  film 
surrounding  a  liquid  oxygen  droplet  by  using  an  assumed  arithmetic 
mean  film  temperature, 


I\  +  T 

_ g_ 


(A-18) 
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APPENDIX  B 

DETERMINATION  OF  A  MAXIMUM  DROPLET  SIZE  FOR  LOX-JET  ATOMIZATION 

The  liquid  oxygen  gaseous  hydrogen  propellnnt  combustion  field  was 
reported  in  Ref.  I  to  have  been  observed  by  simultaneous  motion  and 
streak  photographic  techniques.  Velocity  data  reduced  from  luminous 
traces  on  the  streak  films  exhibited  characteristics  thnt  supported  the 
hypothesis  that  the  highest  velocities  measured  at  any  particular  distance 
from  the  injector  were  representative  of  local  combustion  gas  velocity. 
Considerable  doubt  was  expressed  concerning  the  physical  significance  of 
the  lowest  observed  velocities,  particularly  at  distances  several  inches 
downstream  of  the  photographically  observed  disappearance  of  a  discrete 
liquid  oxygen  jet  spray.  At  approximately  the  end  of  the  visible  liquid 
oxygen  jet,  however,  the  lowest  velocities  approximated  liquid  oxygen 
injection  velocities,  and  further  downstream  more  or  less  abrupt  velocity 
increases  were  noted.  It  was  stated  (Ref.  l)  that  "The  inflections  to 
increasing  lower  boundary  velocities  probably  correspond  with  completion 
of  liquid  oxygen  jet  atomization.  For  a  short  distance  downstream  of  that 
point,  the  lower  boundary  must  reasonably  be  attributed  to  lurge  liquid 
oxygen  droplet  trajectories." 

Droplet  ballistic  equations  show  that  a  large,  massive  droplet  will  be 
accelerated  more  slowly  by  a  given  gaseous  flow  around  the  droplet  than 
will  a  small  droplet.  If  the  lowest  velocity  streaks  correspond  with 
large  droplets,  then  (neglecting  initial  velocity  and  gas  inhomogeneity 
effects)  those  must  be  the  largest  droplets;  otherwise,  even  lower  velocity 
streaks  would  have  been  observed.  This  observation  has  been  used  before 
(Ref.  3)  in  combustion  model  calculations  ns  an  indication  of  the  maximum 
size  droplets  *hat  might  occur  in  a  burning  spray. 
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Even  if  the  lowest  streak  velocities  are  not  known  to  result  from  large 
droplets,  it  can  be  hypothesized  that  if  an  arbitrarily  large  droplet 
were  formed,  it  would  produce  streak  velocities  even  lower  than  those 
actually  observed.  If  that  nrbitrnry  diameter  were  then  successively 
reduced,  a  particular  size  might  be  found  that  would  accelerate  so  that 
its  velocity  history  would  duplicate  (or  approximate)  the  lowest  observed 
velocities.  That  particular  size  could  '■.hen  be  considered  to  be  the 
maximum  droplet  size  produced  in  the  experiment  under  investigation. 

That  concept  was  applied  in  determining  the  maximum  liquid  oxygen  droplet 
size  that  the  atomization  process  model  should  be  permitted  to  produce. 

It  was  assumed  that  the  rocket  combustion  gns  flow  field  could  be  pre¬ 
scribed  throughout  an  appropriate  length  of  the  combustion  chnmber. 
Parameters  of  interest  include  velocity,  pressure,  composition,  molecular 
weight,  density  and  the  transport  properties  (viscosity  and  thermal 
conductivity).  These  were  derived  from  construction  of  an  upper  boundary 
around  the  experimental  streak  velocity  data  and  from  an  assumption  that 
the  combustion  gas  at  every  axial  position  corresponded  to  the  injection 
mixture  ratio. 

A  single  droplet  of  liquid  oxygen  was  then  imagined  to  be  inserted  into 
that  flow  field  at  some  point  and  its  subsequent  history  was  calculated 
under  the  assumption  that  its  evaporation  combustion  had  no  discernible 
effect  on  the  massive  combustion  gas  flow.  Because  this  is  an  initial 
value  problem,  the  calculations  must  begin  from  assumed  values  of  droplet 
diameter,  temperature,  and  velocity  as  well  as  from  a  particular  position 
in  the  chamber.  Here,  interest  lies  in  droplets  of  arbitrary  diameter 
and  having  liquid  oxygen  injection  velocity  and  temperature,  being  first 
exposed  to  the  combustion  gas  stream  at  about  the  point  of  complete 
primary  atomization. 
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The  equations  used  in  this  study  are  the  sane  as  those  used  in  the 
complete  combustion  model,  i.e.  Eq.  (7),  (8a),  (8b)  and  (9)  for  droplet 
ballistics,  Eq.  (10)  and  (ll)  for  droplet  temperature  change,  and 
Eq.  (30)  for  droplet  evaporation.  The  physical  and  transport  properties 
employed  * »re  similarly  evaluated  (as  discussed  in  Appendix  D) .  The 
equations  were  programmed  for  solution  with  an  IBM  7094;  the  solution 
was  based  on  simple  stepwise  progression  down  the  length  of  the  chamber 
from  the  droplet  insertion  point.  Portions  of  the  read-in  data  and  of 
the  computed  results  were  machined  plotted  to  give  direct  visualization 
of  how  closely  the  velocities  matched. 

The  calculated  behaviors  of  six  sizes  of  liquid  oxygen  droplets  with 
initial  temperatures  of  180  R  and  initial  velocities  of  180  ft/sec,  which 
were  assumed  to  have  been  inserted  3*0  inches  downstream  of  the  injector 
into  a  combustion  gas  flow  field  based  on  the  experimental  data  of  run 
417,  PETER,  HIA  (Ref.  l),  are  shown  in  Fig.  B-l  a,  b  and  c.  Comparable 
results  are  shown  in  Fig.  B-2  a,  b  and  c  for  similar  liquid  oxygen  drops 
inserted  4.5  inches  from  the  injector  at  200  ft/sec  into  a  combustion 
gas  flow  field  prescribed  from  run  431,  H2TER,  PRA. 

Examination  of  these  calculated  results  suggests  that  droplets  having 
diameters  between  173  and  200  microns  might  correspond  to  the  maximum 
size  droplets  that  could  hove  produced  the  lowest  velocity  streaks. 
Comparison  with  other  calculated  results,  in  which  droplets  formed  closer 
to  the  injector  were  considered,  and  recognition  that  the  actual  g.ts 
densities  over  the  first  few  inches  of  calculation  must  be  higher  than 
those  used  in  the  calculations,  led  to  the  selection  of  230  microns  as 
the  maximum  liquid  oxygen  droplet  diameter  for  use  in  the  model. 
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APPENDIX  C 

FLAMMABILITY  LIMITS  AND  FLAME  VELOCITIES  FOR  FUEL-RICH 
hyerogen/oxygen  MIXTURES 

The  f lammabi lity  limits  of  hydrogen/oxygen  gas  mixtures  at  low  tempera¬ 
tures  and  high  pressures  and  under  relatively  high-velocity  flow 
conditions  have  not  been  systematically  determined.  Values  for  com¬ 
bustibility  evaluation  in  the  combustion  model  calculations  were  derived 
from  examination  of  the  effects  of  these  variables  on  the  limits  for 
less  severe  conditions  than  those  encountered  in  the  rocket  experiments. 

At  atmospheric  pressure,  values  ranging  from  92  to  96  mole  percent 
hydrogen  have  been  reported  (Ref.  23)  for  the  upper,  fuel-rich  flamma¬ 
bility  limit.  Variability  in  the  limit  composition  arises  from  differences 
in  type  and  size  of  apparatus  employed,  strength  of  the  ignition  source, 
direction  of  flame  propagation,  purity  of  reactants,  and  initial  mixture 
temperature.  A  reasonable  base  value  appears  to  be  9A  percent  hydrogen 
at  291  K  and  1  atmosphere. 

Constantine  (Ref.  24)  has  used  a  concept  attributed  to  Zabetakis  to 
calculate  the  effects  of  renctant  temperature  on  the  upper  flammability 
limit.  The  limit  mixture  is  said  to  occur  at  compositions  such  that 
combustion  product  temperature  remains  a  constant.  The  lowered  initial 
enthalpy  that  results  from  mixing  liquid  oxygen  with  gaseous  hydrogen 
was  calculated  to  narrow  the  upper  limit  (i.e.,  lower  the  hydrogen  con¬ 
centration)  by  only  about  0.2  percent.  Varying  the  initial  hydrogen 
temperature  had  a  much  more  pronounced  effect.  Mixing  liquid  oxygen 
with  hydropen  at  80  K  gave  a  calculated  limit  of  92.8  percent  hydrogen, 
which  broadened  linearly  with  increasing  hydrogen  temperature  to  95.5 
percent  at  580  K. 
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Raising  the  pressure  apparently  broadens  the  upper  limit  for  these 
reactants;  an  increase  from  about  92.5  percent  hydrogen  at  atmospheric 
pressure  to  about  95.5  percent  at  15  atmospheres  has  been  reported 
(Ref.  25).  This  trend  has  been  interpreted  as  supporting  a  thermal  flame 
propagation  mechanism  rather  than  n  chain  branching  mechanism  (Ref.  26). 

The  superposition  of  a  flow  velocity  on  limit  mixtures  of  hydrogen  and 
oxygen  apparently  hns  not  been  studied.  There  is  some  evidence  (from 
comparing  quiescent  tube  nnd  flat-flame  burner  results  for  other  fuels) 
that  suggests  that  broadened  limits  may  be  experienced  (Ref.  26).  For  a 
purely  thermal  propagation  mechanism,  however,  turbulence  of  the  flowing 
stream  may  enhance  local  quenching  nnd  result  in  narrowed  limits  (Ref.  27). 

For  lack  of  more  definitive  information,  values  of  upper  flnmmability 
limit  were  selected  for  model  computations  with  an  assumed  simplification 
that  the  effect  of  pressure  is  offset  by  the  effects  of  flow  nnd  turbulence. 
Thus  values  of  92  to  65  percent  hydrogen  were  used. 

Quantitative  information  regarding  the  turoulent  flnme  velocities  through 
hydrogen  oxygen  mixtures  are  even  less  abundant  than  flnmmability  limit 
data.  Most  investigators  have  measured  laminar  flntne  speeds  nnd  the 
majority  of  the  available  data  are  for  oxidiz^r-rich  or  near  stoicniometric 
mixtures.  For  laminar  flames,  the  highest  velocities  have  been  observed 
with  mixtures  that  are  slightly  fuel-rich,  e.g.,  about  70  to  75  percent 
hydrogen  (Ref.  28  and  29).  These  highest  laminar  velocities  increase  with 
pressure;  data  of  Ref.  28  and  29  for  pressures  of  0.5  to  15  atmospheres 
are  quite  well  fit  by  an  exponential  equation  of  a  form  given  in  Ref.  28, 

S  (cm  sec)  =  1200  +  315  log  [P  (atm)]  (C— 1 ) 
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Flame  speeds  at  90  ntomospheres ,  however,  exceeded  that  given  by  Eq. 

(C-l)  by  about  75  percent  so  that  extrapolation  to  pressures  g  enter 
than  15  atmospheres  is  certainly  queationable .  A  more  often  reported 
relationship  between  laminar  flame  speed  and  pressure  is, 

su  .  a  Pb  (C-2) 

Values  of  the  pressure  exponent  apparently  range  from  about  -1  2  to 

-i-l,  depending  upon  the  specific  renctants  involved  and  on  the  presence 

and  concentration  of  diluents.  Some  investigators  have  suggested  that 

the  exponent  is  related  to  the  order  of  reaction,  n,  through  b  =  1  2  (n  -  2); 

experimental  data  should  then  favor  values  of  b  =  -1  2 ,  0 ,  +1  2 ,  etc. 

The  laminar  burning  velocities  of  Ref.  28  and  29  are  reproduced  in  a 
log-log  plot  in  Fig.  C-l.  It  is  seen  that  an  equation  of  the  form  (C-2) 
can  include  the  data  at  90  atmospheres ,  but  does  not  fit  the  low-pressure 
data  nearly  so  well  as  did  Eq.  (C-l). 

Turbulent  flame  velocities  through  a  stoichiometric  hydrogen  oxygen 
mixture  have  been  determined  at  two  pressures  using  an  extension  of  a 
Bunsen  burner  method  (Ref.  30).  These  two  data  points  are  also  shown  in 
Fig.  C-l.  While  the  turbulent  value  at  14.6  atmospheres  is  1-3  4  times 
the  laminar  value,  considering  them  to  be  ±30  percent  away  from  the 
equation  shown,  which  correlates  the  laminar  dat it  makes  the  difference 
seem  of  dubious  significance. 

In  Ref.  1, however,  the  velocity  data  reduced  from  luminous  streak 
photographs  indicated  that  the  maximum  gas  velocities  just  upstream  of  the 
position  of  abrupt  velocity  increase  were  very  appreciably  higher  than 
the  laminar  or  turbulent  data  in  Fig.  C-l.  This  is  roughly  indicated  by 
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a  circled  triangle  on  that  graph.  Considering  that  symbol  as  a  valid 

turbulent  flame -velocity  point  permits  the  approximation  shown,  i.e., 

1  2 

t’t«  p  .  This  relationship  was  applied  in  the  model  calculation,  but 
the  proportionality  constant  was  selected  arbitrarily  to  force  the  posi¬ 
tion  of  achieving  combustible  conditions  to  match  the  observed  velocity 
jump  position. 

Other  variables  that  are  probably  influential  in  determining  flame  speed 
have  not  been  included,  simply  because  no  data  were  found.  The  laminar 
flame  speed  is  known  to  decrease  sharply  with  increasing  hydrogen 
concentration,  above  about  70  percent  (Ref.  29),  but  apparently  no  deter¬ 
minations  with  near-limit  mixtures  have  been  obtained.  Similarly,  the 
effect  of  lowering  the  reactant  mixture  temperature  is  not  known. 
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APPENDIX  D 

PHYSICAL,  THERMODYNAMIC,  AND  TRANSPORT  PROPERTIES  FOR  THE 

COMBUSTION  REGION 

HIOPERTIES  FOR  EQUILIBRIUM  COMBUSTION 
GASES  AT  VARIOUS  MIXTURE  RATIOS 

The  computed  results  of  frozen  thermodynamic  equilibrium  performance 
calculations  using  the  Rocketdyne  n-element  program  (Ref.  3l)  were  used. 
Properties  were  calculated  for  mixture  ratios  ranging  from  0.3  to  48 
pound  O^/pound  H }  at  330  psia  pressure.  The  properties  needed  in  the 
combustion  model  computations  are  presented  in  Tig.  D-l,D-2,and  D-3  as 
functions  of  gas  mixture  ratio. 

Value*  were  read  from  these  graphs  at  18  selected  mixture  ratios  (10  of 
which  were  between  0  and  9)  and  were 'made  available  to  the  computer  as 
tabular  arrays.  At  times  in  the  model  program  when  gas  properties  at  a 
particular  location  in  the  chamber  were  required,  the  equivalent  gas 
mixture  ratio  at  that  location  was  used  to  enter  and  search  the  table. 
Properties  were  determined  by  linear  interpolation  between  the  tabulated 
values  corresponding  with  the  mixture  ratios  immediately  above  and 
below  the  local  mixture  ratio.  If  the  local  effective  mixture  ratio 
exceeded  the  highest  tabulated  mixture  ratio,  the  properties  at  the  highest 
tabulated  value  were  used. 
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PROPERTIES  OF  GAS  FILM  AROUND  A  DROPLET 

Evaluations  of  the  Reynolds,  Prandtl,  Schmidt,  and  Nusse it* numbers  in 
Eq.  10  and  11  were  based  on  mean  values  of  the  appropriate  properties 
within  the  gas  vapor  film  around  a  droplet.  In  the  m  nburning  region,  the 
properties  were  approximated  simply  by  assuming  the  free-strenm  composition 
and  an  arithmetic  mean  film  temperature.  In  the  combustion  region,  such 
an  approach  is  not  at  all  realistic  because  the  combustion  gases  are 
reactive  with  the  oxidant  droplet  vapors;  much  greater  deviations  from 
free-stream  composition  mny  be  required  before  such  a  mean  temperature 
can  be  effected. 

For  that  reason,  menn  film  properties  were  determined  in  the  combustion 
region  by  the  following  method.  The  vnpors  evolved  from  a  droplet  were 
assumed  to  react  instantaneously  with  the  surrounding  combustion  gases  to 
thermodynamic  equilibrium.  The  oxidant  concentration,  however,  must 
decrease  with  increasing  distance  from  the  droplet  surface.  At  some 
point,  the  composition  can  conceivably  be  such  that  the  local  adiabatic 
flame  temperature  is  the  arithmetic  mean  temperature  between  the  droplet 
and  free-stream  temperatures;  the  mean  film  properties  were  evaluated 
at  that  point.  This  was  accomplished  in  the  combustion  model  program 
by  cross-plotting  some  of  the  desired  properties  against  adiabatic  flame 
temperature  from  the  n-element  performance  program  computations  at  several 
oxidizer-rich  mixture  ratios  and  then  curve-fitting  those  cross  plots. 

The  resultant  equations  were: 

Mf(?  -  (4.4311  x  10-5)  loge  (Tfg)  -  3. 157-  x  10'4  ( D- 1 ) 
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k  .  exp  [(2.979  x  1°",')T.  -  3.5967  [ 


(D-2) 


Oxygen  diffusivity  was  related  to  the  combustion  gas  mixture  ratio: 


fV.cg  ’  (°-0U21  +  °'09  Re*>/Pl  •  “eg 


(0.35444  -  0.007382  R  J/P.  ;  >3.5 

CK  A  c8 


(»-3) 


Specific  heat  of  oxygen  vapor  was  evaluated  by  Eq.  ( A— 7)  using  T^.  . 
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APPENDIX  E 


CROSS-SECTIONAL  AREAS  VS  DISTANCE  FROM  INJECTOR  FOR 
TRANSPARENT  TWO-DIMENSIONAL  MODEL  ROCKET 


The  transparent  two-dimensionni  model  rocket  combustion  chamber,  from 

which  experimental  results  are  used  in  this  report  " . had  a  constant 

1.00-  by  21.63-inch  cross  section  from  the  injector  face  to  12.5-inches 
downstream  of  the  injector,  then  converged  to  a  1.00-  by  13.50-inch 
cross  section  nozzle  throat  at  22.7-inches  downstream  of  the  injector. 

The  initial  chamber  contraction  ratio  was  thus  1.60; . During  the 

nainstnge  portions  of  the  experiments  with  hydrogen  oxygen  propellants, 
however,  melting  of  the  pyrex  walls  began  a  few  milliseconds  after  full 
chamber  pressure  was  achieved.  Small  globules  of  molten  glass  were 
continuously  shed  from  the  chamber  exhaust  so  that  the  nozzle  area  gradu¬ 
ally  increased  throughout  each  test  and  the  contraction  rntio  gradually 
decreased."  (Ref.  l). 


The  effective  contraction  ratio  ef f »  at  the  time  during  a  particular 
firing  that  the  streak  velocity  data  were  reduced  was  inferred  from  the 
apparent  combustion  g;  s  Mach  number  7)1  at  the  beginning  of  nozzle 
convergence.  Isentropic  acceleration  of  a  completely  reacted,  frozen- 
equilibrium  combustion  gas  stream  from  a  cross-sectional  area  A^  to 
sonic  velocity  at  the  nozzle  throat,  with  area  A  ,  is  given  (Ref.  32)  by, 
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*Aeff  .A 


M 
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The  melting  and  erosion  of  the  Pyrex  was  seen  photographies  ny  to 
occur  fairly  generally  downstream  of  about  7  inches  from  the  injector, 
both  the  motion  pictures  and  posttest  inspection  suggested  that  the  major 
configuration  changes  probably  occurred  in  the  nozzle  section.  The  chamber 
cross-sectional  area  per  injection  element  was  calculated  for  combustion 
model  input,  therefore,  with  a  wall  erosion  factor  applied  only  in  the 
nozzle.  This  factor  was  assumed  to  have  a  value  of  unity  near  the 
beginning  of  nozzle  convergence  and  to  increase  linearly  to  a  maximum 
value  of  1.60/«^^  at  the  nozzle  throat. 

The  equations  u  ed  were: 

A  =  A  .  ,  0  <  x  *  12.50  (E-2) 

CC  CCl  '  7 

A  =  F  A  To. 538  +  0.0924  J  25  -  (x-^.SO)^1  ]  , 

cc  erosion  cci  L  V  '  7  J 

12.50  <  x  <  14.82  (E-3) 


F  A 

erosion  cci 


[ 


1.059  -  0.0483 


(x 


12.50)]  , 


14.82  <  x  <  20.47 


(E-4) 


A  =  F  A  1.741 

cc  erosion  cci  L 


0. 1482  v  25  -(22.7-x)^  j 


> 


20.47  <  x  <  22.70 


where: 

(1-6  -•*,,«>  (*  -  12-50) 

’erosion  °  1  *  ..  (22.70  -  12.50) 


(E-5) 

(£-*>) 
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NOMENCLATURE 


a 

A 

A 


eff 


B. 


CA 

S 

Ct 
D 
D 


20 


D 


30 


D_ 


32 

r 

V 

e 


*div 


*1 

F(n) 

«c 

II 


velocity  of  sound,  ft,  sec  or  an  arbitrary  constant 

area,  sq  in. 

surface  area,  sq  ft 

effective  chamber  contraction  ratio 

spray  characterization  parameter  or  an  arbitrary  constant 
spray  orop  size  parameter  (Eq.  l) 
arbitrary  constants 

specific  heat  at  constant  pressure,  Btu  lbm-R 

spray  atomization  rate  parameter 

drag  coefficient 

Craya-Curtet  numb  *r 

diameter,  inches  or  microns 

surface-number  mean  droplet  diameter  for  a  spray 
volume-number  mean  droplet  diameter  for  a  spray 
volume-surface  (Sauter)  mean  droplet  diameter  for  a  spray 
diffusivity  of  one  gas  into  another,  sq  ft,  sec 
base  of  Naperian  logarithms  =  2.71828 
functional  symbol 

fraction  of  spray  element  gas  flow  diverted  out  of  the  spray  combustion 
zone 

fraction  of  initial  gas  momentum  still  associated  with  gaseous 
propellants 

factor  for  suppressing  atomization  j”«t  downstream  of  the  injector 
gravitational  constant  =  32.17  lbm-ft  lbf-sq  sec 
heat  of  vaporization,  Btu/lbm 
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It 

I  =  momentum,  lbf 

k  =  thermal  conductivity,  Btu/ft-hr-R 

K  =  a  defined  droplet  drag  parameter  (Eq.  7) 

=  chamber  length  to  accomplish  turbulent  mixing,  inches 
M  =  molecular  weight,  lbm  lb-mole 

=  mach  number 

' 

N  =  number  flowrate  of  spray  droplets,  sec 

Nu  =  Nusselt  number 

Nu^  =  2  +  0.6  Re1  2PrX  3 

Ni^  =  2  +  0.6  Re1  V'  3 

p  =  spray  mean  drop  diameter  variable  (Eq.  5**) 

P  =  pressure,  lbf  sq  in. 

Pv  =  liquid  vapor  pressure,  lbf/ sq  in. 

Pr  =  Prandtl  number  =  3600  c  fj^k 

q  =  spray  mean  drop  diameter  variable  (Eq.  54) 

r  =  radius,  inches 

R  =  universal  gas  constant  =  1544  ft-lbf  lb-mole-R  or  gas  mixture  ratio, 

if  subscripted 

Re  =  Reynolds  number,  p  D  U/(l2  \i) 

Sc  =  Schmidt  number  =  fi/(p  £) 

S^  =  laminar  flame  speed,  ft  sec 

t  =  time,  seconds 

T  =  temperature,  R 

u  =  an  integrating  parameter  (Eq.  59  ) 

U  =  velocity,  ft  sec 

U  =  kinematic  or  space-mean  velocity,  ft  sec 

K 

v  =  cumulative  volume  fraction  of  spray  smaller  than  diameter  D 


w 


droplet  weight,  lbm 
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=  weight  flowrate,  lbm  sec 
=  local  spray  production  rate,  lbm  sec 
=  distance  along  the  direction  of  flow,  inches 
=  incremental  distance,  inches 
*  mole  fraction 

=  a  dimensionless  droplet  diameter  function  (Eq.  58) 

=  heat  transfer  factor  (Eq.  11 ) 


Greek  Letters 


y 

6 

t r 

P 

o 

r(P) 
r  (p) 

Uj  r 

L 


ratio  of  specific  heats  (c^  c^) 
spray  characterization  parameter 
viscosity,  lbm  ft-sec 

3.14159 

density,  lbm/cu  ft 

surface  tension,  lbf  ft 

Ganna  function  =  J°°  u^-^  e  U  du,  p  £  -  1 

incomplete  Gamma  function  =  J  u^-^  e  U  du,  p  ^  -  1 
sunnation  symbol 


Superscripts 


(pri 


ime 


\=  gas  properties  evaluated  at  an  effective  atomization  distance  away 
from  the  liquid  jet  surface 
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* 
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Subscripts 


« 


A 

b 

cc 

eg 

cr 

d 

def 


est 

f 

final 

K 

gs 


i 


j 

l 


max 

mean 

n 

o 


0 

P 


2 


q 


local  surface  atomization 
gas ,  spray  zone  boundary  value 
combustion  chamber  value 
combustion  product  gas  value 
critical  state  value 
droplet  value 

value  within  the  mixing  zone  between  a  gas  jet  and  its  surroundings, 
referred  to  as  a  velocity  defect  zone 

estimated  value 

value  within  the  gas  film  surrounding  a  droplet 
value  at  completion  of  combustion 
value  associated  with  a  gas 

value  associated  with  the  non-sp~ay-bearing  gas  stream  in  the 
combustion  region 

gaseous  hydrogen  value 

initial  value 

jet  value 

liquid  (subcritical  temperature,  condensed  state)  value 
maximum  value 
average  value 

value  associated  with  the  nth  Bpray  group 
stagnation  value 
oxygen  value 

value  at  or  preceding  the  position  furthest  from  the  injector  to  which 
a  portion  of  a  gas  jet  penetrates  with  injection  velocity  maintained 
or  spray  mean  drop  diamete  variable 

spray  mean  drop  diameter  variable 
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V 

X 

1 
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1-2 


»  relative  value  or  value  at  an  arbitrary  radius 
=  spray  zone  value 
=  turbulent  flame  value 
=  vapor  state  value 
=  wall  value 

=  value  at  an  arbitrary  axial  position 
=  value  at  the  beginning  of  n  A  x-  increment 
s  value  at  the  end  of  a  A  x-  increment 
s  value  produced  within  a  A  x-  increment 
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TABLE  1 


INITIAL  MODEL  CONDITIONS  FOR  CASES 
COMPARABLE  WITH  EXPQUMINTAL  DATA 


Initial  Value  of 

Experiment  for  Comparison  (Ref.  1) 

Parameter 

Test  No.  417 

Test  No.  431 

Pj,  psia 

457 

455 

A  . ,  sq  in.  element 
cci  ' 

0.483 

0.525 

1.312 

1.344 

. ,  lb  /sec-element 
02i’  m 

0.592 

0.655 

\i’  " 

185 

185 

AQ  aq  in  ./element 

0.01368 

0.01368 

U  ft/sec 

2 

89.6 

99.1 

w„  lb  /sec-element 

l^i  nr 

0.1200 

0.0921 

VR 

~  465 

~  305 

A^^,  sq  in. /element 

0.03452 

0.03452 

U„  ft  sec 

Hr,* 

2670 

1360 

R. 

l 

4.93 

7.11 
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TABLE  2 

StMMAHY  OF  EXPERIMENTAL  CONDITIONS 
AND  VISIBLE  LOX  JET  LENGTHS 


“ 

Test 

No. 

Injector  End 
Pc,  psia 

Mixture 
Ratio , 
o/f 

Hydrogen 

Injection 

Temperature, 

R* 

Approx.  Visible 
LOX-Jet  Length, 
inches 

417 

[470]** 

4.9 

-  500 

2.8  to  3.2 

418 

516 

5.0 

[chilled'1 

3.5  to  3.8 

419 

516 

5.0 

[chilled] 

3.5  to  3.7 

431 

[430]** 

7.1 

315 

4.1  to  4.3 

432 

516 

4.1 

205 

3.0  to  3.2  ! 

433 

522 

5.2 

210 

4.1  to  4.3 

432 

459 

5.8 

200 

3.8  to  4.1 

*The  temperatures  measured  in  the  injection  manifold  were  con¬ 
sidered  in  the  combustion  model  to  be  stagnation  values. 

**  Extrapolated  from  downstream  pressure  data  by  comparison  with 
measured  values  from  the  other  tests. 
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Figure  1.  Schematic  Representation  of  the  Combustion  Processes  for  Coaxial  Jet  Injection  of 
Gaseous  Hydrogen  and  Liquid  Oxygen  (Ref  l) 
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Figure  4.  Calculated  Nonburning  Region  Structure  for  Conditions  of  Test  417 
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Figure  5.  Coaxial  Jet  Injection  of  Liquid  Nitrogen  Into  Ambient  Air  (Ref.  14) 


Figure  6 


Photograph  of  Water  Jet  Atomization  From  Two 
Coaxial  Jet  Elements  W'th  Gaseous  Nitrogen 
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BEING  PRODUCED  SPRAY  BEING  PRODUCED 


Figure  8.  Weight  Distribution  of  Spray  Being  nroduced 


Figure  9*  LOX— Jet  Atoaization  and  Spray  Vaporization  in  the  Nonburning  Region 


LOX  JET 
SURFACE 


(b)  TEST  NO.  431 


Figure  11. 
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Oxygen  Concentration  Profiles 
in  the  Nonburning  Spray  Zone 
at  Various  Distances  From  the 
Injector 
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Velocities  of  the  Spray  Zone 
Boundary  Gases  and  Selected 
Droplet  Size  Groups  in  the 
Nonburning  Spray  Zone 
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Fiq^ire  13.  Spray  Element  Droplet  and  Gaa  Temperatures 
in  the  Nonburning  Region  (Test  417) 
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DROPLETS  GASES 


Figure  14.  Maximum  Equilibrium  Oxygen  Concentration  in  Hydrogen— Oxygen 
Gas  Mixtures  That  Contain  LOX  Spray 
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Figure  15.  Calculated  Craya-Curtet  Number 
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Figure  16.  Calculated  Combustion  Region  Characteristics,  Test  417 
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Figure  17*  Comparison  Between  Experimental 
Data  and  Calculated  Velocities 
and  Pressures  of  the  Best-Fit 
Cases 
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DISTANCE  FROM  INJECTOR,  INCHES 


Figure  18.  Progress  of  Oxygen  Atomization  and 
Reaction  in  the  Combustion  Region 
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Figure  B-l.  Calculated  Behaviors  of  Variot 
Sizes  of  Liquid  Oxygen  Drople 
in  a  Combustion  Gas  Flow  Fiel 
Prescribed  from  Bun  *il7,  PETE 
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Figure  B-2.  Calculated  Behaviors  of  Various 
Sizes  of  Liquid  Oxygen  Droplets 
in  a  Combustion  Gas  Flow  Field 
Prescribed  from  Run  431,  PETER, 
PRA 
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Figure  D-2»  Equilibrium  Molecular  Weight  and  Hatio  of  Specific  lleata  for  1  jOX/lU 
Combustion  Product  Gases  at  530  Psia 
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Simultaneous  equations  describing  rocket  propellant  injection,  atomization, 
mixing,  vaporization,  and  combustion  are  formulated  for  a  cylindrical  liquid 
oxygen  jet  surrounded  by  an  annular  gaseous  hydrogen  stream.  The  formulation 
is  basod  on  division  of  the  combustion  chamber  into  a  nonbuming  region  neart 
the  injector  and  a  combustion  region  further  downstream.  Processor  in  the  non- 
burning  region  are  required  to  satisfy  certain  combustibility  criteria  before 
combustion  is  possible.  The  propellants  then  pass  through  a  plane  flame  front 
into  the  combustion  region.  A  simplified  treatment  of  turbulent  aixihg  between 
unlike  reacting  gas  streams,  one  of  which  contains  a  liquid  oxygen  spray, 
characterizes  the  combustion  region  model.  Solutions  of  the  system  of  equations, 
obtained  wiwh  a  digital  computer  program,  are  detailed,  discussed,  and  compared 
favorably  with  experimental  information  froa  Part  I. 


DD  .MR.  1473 


racussimo 


